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FORE VeD 


Thig technical report is preparad for pretentation at 


training seminars modeling plate Glementa in WASTREAN CHASA 


STRuctural ANalysis Program). Thiga training weminar is weausally 
a © given at NASTRAN User Colloquiums or NASTRAN Applications Short 


Courses. 
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NASTKAN (NASA STRuctural gWalysia Program) i the mogt 
widely used gtneral purpose atructural analysis program in the 
world. This program wax originally developed in the mid-gixties, 
and ¢he firat version was relwaged in 10@8., Th program 
davealuped was gsponored by the National Aercuautics and Space 
Administration (NASA). This development in the form of 
enhancements and maintenance was continued until] 1072 as a wsing.a 
Government version under the sponsorship of NASA. Sicce then the 
Government vergion ia being deveioped as COSMIC-NASTRAN, while a 
commercial version is being marketed by McNea! and Swendler 
Corporation (MSC) aa MSC-NASTRAN. Even though the basic 
atructure of these two versions remaina the game, there are 
waignificant differences in capabilities and efficiencies. Around 
1989 the MSC introduced a new plate bending alement called the 
QUaDd to its element library. {t igs one of the most veraatile 
plate bending elements, even though its theoretical basie ia 
Somewhat controversial. The basic atrength of this element is 
thal i¢ was wobigcted to excluaive numerical testing, and a 
number of empirteoe)] aviuatmenta wo re made to conform the regulis 
to Brown #oftbiona. jhe QUADS @#lement embodias a number of 


Ampoovements over lat @arlier elementa (QDMUN) ODMEMI, QDMEMZ. , 


:: It is an isoparametric formulation. 

b. ft models ftnplane (membrane) behavior more realistically. 

c. The layered composites modeling capability ig extensive. 

ad. Membrane-bending coupling can be modeled pealiatically. 

e. It i# the only plate element with an offset feature. 

> Sipe It ig a convenient element for modeling laminate plates. 

@. The same element can be used in modeling sandwich plates; 
even though it is not as simpie to model sandwich plates 
with composite face sheets. 


h. A single piate element replaced ail other elements. 


Until around 1987, COSMIC~NASTRAN did not have a similar 

(QUAD4) capability. The absence of this capability rupre#zentead a 

‘gnificant inconvenien:se, in particular, for uwodeling layered 
cen pogzites. Betwees 1983-1966 a QUADS element was dovelopsd 
under the uponsorahip of the Air Force (Flight Dynamics 
Laboratory) for use in the program ASTROS (Automated $tructural 
Uptimiecation Evetem). At the w-me time this slement was 
incorporated into COSMIC~-NASTREAN The ASTROS-QUAD4 element if 
wimilar to the MSC-QUAD4, but there are significant differences 
in the theoretical formulation between the two elements. 

This report provites an informal background for modeling with 
the QUAD4 @lement. Ut contains over 15 problems to illuatrate 
Various options of the wlament However, the report was intended 
only as background material to be uses in conjunction with a one~ 
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APPENDIX A 


THE QUADS ELEMENT 


This appendix provides the theoretical development for the QUAD4 
element that has been installed into ASTROS. An overview of this element is 
given in Subsection 5.3.3, while this appendix provides detailed information 
on the element. This detail is necessary because, unlike the other elements, 


the ASTROS QUAD4 element has not been documented elsewhere. 
A.l DISPLACEMENT FUNCTIONS 


The QUADS element inas two distinct element coordinate systems. These 
are the "user defined" element coordinate system as defined by the element 
connectivity data and the “internal element" coordinate system, which is 
defined as having its origin at Gy (X°¢, Ys, Z°p). This origin is computed 
by taking the average of the grid point coordinates. The positive X- and Y- 
axes of the internal element coordinate system are defined with the aid of two 
points, Gyp and Gyp described below. 


Vi3 and Vou are defined as the unit diagonal vectors as illustrated in 
Figure A-1. Thus, the coordinates of points Gyp and Gyp are given by the 
following: 


Cy ow X° + x’ » (Yo +¥), (29 
E {( : E ( E E ( 2) 


é (A-1) 
= ((X°- ¥'), (Yo 4# x"), (22) 
YE (¢ : ~ ( E E ( E’ 


where, X°s, Y°¢ and 2°, are the coordinates of the origins of the internal 
coordinate system and X'p and Y's are the components of the bisector vector of 


the unit diagonals V13 and V2,z. 


The coordinates of points Gp, Gyp and Gyg, are used to define the 
transformation from the internal element coordinate system to the coordinate 
system in which the grid points are defined. The internal element coordinate 
system is necessary to correctly handle irregular-shaped and non-planar 


elements and is henceforth referred to as the “element” (E) coordinate system. 
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Using 2-D interpolation functions, the geometry field at any point (€, 
n) in the element cross-section (see Figure A-2) is defined, where the nodal 
curvilinear coordinates are related to the nodal cartesian coordinates system 
in the element coordinate system by the following relationship: 
4 i 
(XE(Em)) = jE, Nan) (XB) 


where i refers to grid point i, and 


1 
{X_) = (Xp, Ye, Ze) at node i, 


N(é,m)) are the interpolation (shape) functions which define the con- 
tribution of each node at a given point with the element. 


These functions and their derivatives are: 


Ny om 1/4(1 + €€4)(1 + nq) 

aN 

—i -« des + nz) (A-2) 
én 4 


The deformations of the element are also represented with the identi- 
cal interpolation functions: 


4 i 
{Up(€,7)} = rom Ny (€,7) (Ug) (A-3) 


where (utp) - (Up, Ve, We, Pxe, *yE, Pzpj)T represents the vector of 
displacements at grid point 1 in the element coordinate system. 


A.2 STRAIN-DISPLACEMENT RELATIONSHIP 


The QUADS element incorporates a reduced solid theory for thick 
shells. According to this theory, the element has five dof at each grid, _ 
defined in a coordinate system whose X-Y plane is tangent to the mid-surface 
of the shell at the given grid point. The z-axis, therefore, is the normal to 
mid-surface at that point. In our nomenclature, this is called the "C* systen 
(Figure A-2 and A-3). | 


A generalization of the "C" system, called "I" system, incorporates 
the characteristics of the "C" system at a general point on the mid-surface of 
the shell element, normally the integration point (Figure A-2). 
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In order to establish a common definition for "I" and “C” systems, 


consider the following steps: 


(A) The tangents to mid-surface at a given point (€,”) are: 


(B) 


(C) 


| 

8 yy 

z 4 an, [* i 

(1) = ——-E - = ity (A-4) 
ae i-1 8€ |2J, 

8 4y 

z 4 8N x i ; 

(Yeo) = ———E - gc tty (A-5) 
én fel 8 |p E 


i 
x 
(| are the coordinates of grid points in “E" system. 
ZJE 


The axes of the new system then follow: 
(Ve } x (Ve) 
pone SERS 


(2}y7¢ = (Vp) = 
[(Ve,) * (Vey) 


(Y). x (2) 
(yo o —E____1/e_ (A-6) 
je x (Z} r/c] 


Wyo = (2)ryy¢ * (Xr/c 


Finally: 

[TIE] = [{X)q{¥)q(Z)z)T (A-7) 

{TcE}i = pcxyleyytezj4jT (A-8) 
¢c c ¢ 


Note that the "C" system is not necessarily invariant when we go 
from one grid to the next. This is due to the possible warping 


of the element. 


Since the ultimate goal of this discussion is to establish a relation- 


ship between the element strains (which are defined in the "I" system), and 


the nodal displacements (defined in the "E" system), it is necessary to 


develop a series of transformations along with the strain-displacement rela- 


tionships. 


3c 


Consider the five dof's in the "C" system at each grid point "i" to be 
arranged in the following manner (Figure A-3): 


i uj} i a)! 
(Uy = Av > (A-9) 
c a c 
fed c 


In order to be compatible with the other dof’s in the model, these are 
related to the six dof at that grid point, defined in the *E” system, by the 
relationship: 


i i 
(U) = [tTcEt(u) 
c E 


i 010 4 
(@) = {TcE]1(6) 
c -100 E 


(A-10) 


The extra transformation in the rotational case is a result of the 
difference in the definition of rotations for "C" and "E" systems (Figures A-3 
and A-4). 


The same five dof’s are related to six dof’s in the "I" system by 
using the transformations developed in Equations A-7 and A-8. Considering 
Equation A-10 and A-3; 

(U) [TIE] : N terce)$) "(uy [TIE] : N (vy 
i-1 ' c inl} E 


4 i 
= = Ny(T)](U) A-11 
i-l i(T] E ( ) 
and 


4 T i 4 T 
(8)y = [TIE] ENy{(TCE)!) (6) = [TIE] E.Ny[(TCE)4L) 
i=] c {el 


(A-12) 


010 { 4 1 
[TCE}4(e)) =  S Ny{a}t(e) 
-100 E- i=l E 


Note that while [T] is invariant, [A] depends on the direction of the normal 
to mid-surface at each grid point. 
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At a point along the Z-axis of "I" system, at a level of 2 = (ty79, 
where, 
4 
t - = Nyt 
I Poor 


is the thickness of the element evaluated at this particular integration 
point, the dof’s in "I" system may be written in the following form: 
(uyly = (uly ; (ugly = $tz/2tely (A-13) . 
The strain-displacement relationships can now be developed, using 


these rearranged dof's: 


tx )M sufax afax 0 0 

Ce OS ey av/ay - | 0 @/ay Ol(uy)y — (A-14) 
IxyJI du/dyt+dv/dx) x a/ay &/ax 0 
€x ]B gu/ax 8 a/ax 0 0 

(ep)r ee ia . év/ay - |} 0 @/dy O}{ug)y  (A-15) 
IxyJ I du/dy+8v/dx) 7 a/ay 8/ax 0 


yz aw/ay+dv/az 00 a@/ay| 0 8/az 0] (uy 
(ys) = - ~ | (A-16) 
Izy 1 aw/8xtdu/dz 00 d/adxja/dz 0 Of lupsy 
Inserting Equations A-11 through A-13 into Equations A-14 through A- 


16, and considering the following: 


2 (up)y - 2 zl6)y = (Hy 
oz 8z 


and (A-17) 
a/ax aN; /dx 
&/ay ~- 5, 8N; /dy 
1 Ny 


we arrive at the following general relationships: 
8N;/ax 0 0 
i 
{em}y = £ 0 6N;/dy 0 Ee (A-18) 


im} 
aNy/ay 8Nz/ax 0 
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6N;/ax 0 0 


cr, 4 . 4d 
(ep=)y = on 0 — 8Ny/ay O | [A] (8) (A-19) 


8N;/8y aNz/ax 0 
4 [00 8N;/ay } 0 Ny O)]f {[T] | 0 vu yi 
(ws}1 - {2 1 sides lapdaher oon (A-20) 
f-1 [00 aNy/ax | Nz O OF! O | [alt ale Je 


or, collectively: 


ee 0 oO | u yi 
| | 
€M O 8N;/8y O| 0 os 0 Vv 
| 
8N;/dy 8N;/ax O| | w 
eo lew ce cote es ceed ate | Silee eid goewteeeeeie eels tems J--- 
. |} [angsax 0 O ! e 
4 let | ‘ 
(elt = yeBr =, 0 i O aNy/ay 0} |] 0 [LA]*] Jey 
| | 
| ee 8N;/8x 0 O2JE 
ee es Se | w meee cceccncaccccnce 
0 O <aNy/ay | O Ni 0 
%s | (A-21) 
I 0 O @Ny/dx | Ny 0 0 


Since the shape functions Ny are defined in terms of the curvilinear 
coordinates (£€,n), the shape function derivatives are related to the corre- 
sponding Cartesian derivatives in the element [E] coordinate system, by using 


the rules of partial differentiation, as: 


aN; /8€ 8x/ak Sy/a€ 8z/aé ]{ aNz/ax 
aNy/an } = | 8x/8n By/an 8z/an |4 aNz/ay (A-22) 
any /a¢ 8x/as By/ag az/as || aNy/az 


The first and second rows of the transformation matrix (or Jacobian 
matrix [J]) are the tangent vectors to the surface r = constant and the third 
row is the interpolated values of the nodal normals. (Note the nodal normals 
are evaluated by carrying out the cross product of the two tangent vectors at 


the node point.) 


From Equation A-7 the coordinates in the "I" system are related to the 
coordinates in the "E" system by the following: 


(U)y = [TIE}(U)g 


35 


Yherefore, the derivatives are given by: 
ONy /dx Nf 2€ 
Ani /dy fo ™ 16) 9 ONy/8y 
ON, /dz aN, /8s 
where (A~23) 
#11 $12 0 


(¢] = [TE] (J]"4 = | go, g22 © 
$313 $32 $33 


Note that 8N;/d¢ and 4Nq/)z will be sero when the interpolated normal at the 
integration point coincides with the normal to the mid-surface; e.g., in the 
case of the flat plate (¢3) and $39 are zero). The zero terms in [6], f.e., 


613 and $93. reault from dut products of perpendicular vectors. 
A.3 STRESS-STRAIN RELATIONSHIPS 


Stresses are velated to the previously defined strains by the elastic- 


ity matrix [G}] (where [G) 4s partitioned to give gemuvate mewhrane ectresser). 


oM [ G, {0 40 ey} foams 4 
i dati ee eed eee ae ce eee 
oR a ee es a Sf S87 
Pere pr mee aes fen Jae | fee ratte! 
L res GC re {63 5 Urs Imes e I 
or CA- 26) 
fo}y om (G}yCCesnec ~ Cepdy 
where 
lop) Menbrane stress vecter 
(op) Bending stress vector 
trro) Trensverse ghear stress vector 
{Gq] Membrane sedulf matrix 
Gg] Bending wodult wstrix 
[Gay Teansvetee sheay woduil metiwa 
and gabser (prs SREC* and "TP refer to mechanica! and therwel, veapect ively. 
Sti 


2 ‘sind MRSA SEU AEN eva reece eae lb Bho tae 
cumtneeniatt {ememne =USE A SEN MA NER A SLSOREAUN. eit 


The membrane-bending coupling moduli matrix [G4] will be incorporated 
into the [G] matrix following this discussion of the uncoupled matrices. 


All anisotropic, orthotropic and isotropic material properties sare 
supported. The elastic modulus matrix [GJ is defined in the material coordi- 
nate system and transformed into the user defined element coordinate system by 
means of a transformation angle, #4, which references the user defined element 
X-AXIS or the material coordinate system ID (MCSID) specified by the user. é@y 
is in the X-Y plane of the element as shown in Figure A-5. 


The elastic modulus matrix in the element coordinate system is: 
[Gly = (U}T{6)yIU}] (A-25) 


(Note that since the projection of Xz onto the Xp-Yg plane is parallel to Xf, 


no extra transformations are required between the "E” and "I" systems.) 
The transformation matrix for [G1], [G2] and [Gz] is: 
cos7 by sin? éy coséysingy 
[Uz] .- sin26y cos*éy -coséysingy (A-26) 
—2sinéycosiy 2siniycoséy cos7éy-sin2)y 
and the transformation matrix for [G3] is: 
| coséy sing | 
[U2] = (A-27) 
—sinéy coséy 
For i{sotropic materials: 


(A) Membrane 


{G)] = a) vy 0 (A-28) 
@-y -y 
SYM iy 


(B) Bending 


4 
fe » t. [6 % 
[Co] ET [G3] (A-29) 


a7 


{C) Transverse Shear 


yf dot. re) 
pe 4 
{Gy] ~ i (A-30) 
~ 2K 


where E fs the Young’s modulus; t is the element thickness at the correspond- 
ing integration peint, » is the Poisson's ratio and s/t is the transverse 


shear factor. 


Note that in matrix [G3], the factor "K" is intreduced to compensate 
for the difference in shear distribution though the thickness, which is 
parabolic and not constant as indicated by the displacement function. ‘the 
value of Kel.2 is the ratio of the relevant strain energies. The A; factors, 
which are derived mumerically, are introduced to compensate for the "locking" 


of the element due to excessive shear stiffness. 
For anisotropic materials: 
(A) Membrane 
G6) 2 G23 
[Gj = G22 623 (A-31) 

SYM 633 

{(B) Bending 
% 


Cay ee ok (A-32) 
[Gj oat [G1] 
(CC) Yransverse Shear 


Cy, Gy | 
(Ae? 


wd 
iat 
~~ 


t 
[G3] = ‘| 
© £ Gz Go J 


For orthotropic materials: 


CA\ Membrane 


ae Sarl 


PE, #1 2F2 G 
Levyovoy 


+ 
{ 
Eo 0 : Ca 34) 


Sv G 12 ¢ t. My ava] 4 


2 ttm 4: ashen ht A debe BA SS AAR EU 


aN A: AR OMS RIAN OS SUI ra be een 


A SLATE GA: RPT A Pa SoD, oH wae ee wed Crs ee ett] 


(Bb) Sending 


3 
Go| = kK [6 -3 
[G2 Ti {Gj} (A-35) 


(GC) Transverse shear 


Giz «0 
(A-36) 


t 
[63] we ht 
t 0 Go 


where Ey and Eg ar: the Young's moduli in the principal material axes, v1 {s 
the major Poisson’s ration; Gj2 is the in-plane shear modulus, Gj, and Go, are 
the out-of-plane shear moduli and *s/t 4s the transverse shear factor. 

The derivation of the [G,] membrane-bending coupling matrix begins by 


denoting the strains at the mid-surface as: 


o 
x 

(¢) o ; 
(*m)} ‘a *y (A- 37) 
6 


and the out of plane curvatures as: 


| Kx 
(K) oo | Ry (4-38) 


| Kyy 


Therefeve, the etrains at e@ distance z above the mid-surface of the element 


are: 
6 
ve} om few) - 2{K} (&-39) 
The corresponding 2-D stresses are: 
© 
to} = (G]1( fey) - £UK)) (A-40) 
where [G}y is a (3:3) matrix of elestic mocull. 
The ferces and moments per unit Length are therefore given by: 
th ae 
(Flo @ fF (eo) dz = fF [G}y(te®) - x(k) az (A-41) 
Ze 24 
(F} erey eek 4.07 (G7 


Xb Zp 
(M} = fo (oe) zdz = ff [e}y(-ze%) + 224K) a2 (A-42) 
Ze Za 


(My om e2(Gg pie?) + Lica] tk) 


where t is the plate thickness and I ise the bending inertia. Assuming a 


linear variation of clastic properties between top and bottom surface. 


79 6 +6 
[Gx] = 4 ff © Gaz o oo 8 (A-43) 
: t -t/9 2 
t/ 3 
[Go} = 2 ff ~ Gaz = (£~y[6}} (A-44) 
t/o G.-G 
[Gg] «- 4. f (-z)Gdz = - |B (A-43) 
t? -t/9 12 


Note that the membrane-bending stiffness coupling terms vanish for a eiement 
whose elastic properties are symmetric relative to the wean plane of the 


element. 


By assuming that the elastic modulus has a linear variation between 


the top and bottom surfaces, define: 
G = G) +6/2(Gy - Gg) (A-46) 
Therefore, from Equations A-31 and A-32: 
(A) Membrane 
Go = Gy) + §/2(-126,) (A-47) 
Go= Gy - 66, (A-48) 


(B) Bending 


3 
Gt 
G) = te ee 
121 
(A-49) 
G ~ Zl GF - 656, 


Ma rix [Gj] fe not affected since transverse shears are assumed to have no 


coupling action, 


Therefore, the stress-strain relationship, sliowing for membrane, 


bending, transverse shear and membrane-bending coupling is: 


oN Gy 1 Gy-6¢G, | © €M 
teevosl) ‘iPedowocccullecesocacs |---- ee 
eTor ¢ = | S1-6G, | G2 | 9 B (A-50) 
decectap hecwome son es eer ease 
TS oO 1 O | § &3 TTS 
where 
Cy 
(oy) = ey » Membrane stresses 
xy 
{ Oy 
lepoT) = i ey » Total membrare and bending stresses 
Yxy JTOT 
Tyz 
(ors) = ; Trausverse shear stresses 
Txz 
( 
(ey) «= } ty » Membrane strain 
\ Ixy 
ex 
{ep) = €\- : Bending strains 
“ny JB 
Vyz 
tyts) = , Transverse shear strains 
xz 4 
A.G STIFFNESS MATRIX 


The element stiffness matrix is derived by minimizing the total 
potential energy and is given in sumerical form by mwploying the Gauss-quadra- 
ture integration wethod: 

[K}_ = EEE [BIT [C} [BWW Uedet [5] (A-51) 
where (€, , ¢) are the Gaussian integration point coordinat:s and We, Wis and 
We are the associated weight factors. Det [J] represents the whysical volume 
of tne element as calev‘tated at this point, B 4s the stiati dlaplacement 
yeletionship of Equation A-?] and G is the atress etrain relationship ef 


Equation A-S0. 


Each element a.:iffness matrix partition in the element coordinate 


system, [KyjJgg, is transformed to the global coordinate system by the follow- 


ing transformation: 


T 
[Kiglg o (TEC): (Kgs) peiTEG), (A-52) 


where [TEG], is determined by relating the element cocrdinate system to the 
giobal coordinate system fer grid 4 through the basic coordinate system: 


[TEC]; = [TEB],{TBC}, (A-53) 
A.5 CONSISTENT AND LUMPED MASS MATRICES 


The consistent mass matrix terms are evaluated, neglecting the rota- 
tional inertias associated with the a and § degrees of freedom, by the follow- 


ing expression: 
4 
Mis = E NGNy plJ[ty (A-54) 
n=] 


ie - where Ny ds the shape function for node 1, p is the mass per unit volume, 13} 
{1s the physical area of the element and t, is the element thickness at the 


integration point. 


The lumped wass matrix, which is calculated at the pseudo center 
flle., the average of the element grid coordinates), is prorated to the edges 


based on the distance of the pseudo center from each edge. 
The terms oxy the lumped mass matrix sre evaluated using: 
4 
Miy > 2 Ny pl3[t, (A-55) 


The transformation of the tass matcix to the global coordinate system 
is carried out using the sare transformation matrices as used for the stiff- 


ness matrix in Equation aA-52. 
A.6 STRESS RECOVERY 


The element stresses In partitioned form from Equetion A-50 are 


ts 


orca ramen ata sanders tabetha mares 
SS ORNs Anda lc RS > Tam 
nC ee ee 


weoewevceti fLwesse2xeseef{eeaerereerererirreer see 


TTS 0 G3 ysJmec Ts 
or 
oy G, | G,'} 0 eM en 
Si) Sethe Odea he: aan ase) 
TTS 0 | 0 G3 rrsJMEC 0 UrtsJT 


For a specified grid point temperature, the thermal strain vector is: 


€x 
(emir = ty = {az}(Ty - To) (A-57) 
T 


where {az} = [U}-2 (ay) is a vector of thermal expansion coefficients in the 
element coordinate system. [U] is the strain transformation matrix given in 
Equation A-26 and (ay) is the vector of thermal expansion coefficients in the 
material axes. Ty and Tp are the specified grid point temperature and mid- 


surface (stress-free) temperature, respectively. 
For a thermal gradient T’, the thermal strain vector (€B)y is: 


(Bly = tag (& 7) (A-58) 
For thermal moments (M)t, the thermal strain vector (‘B)ry is: 
(Bigs Hye tGaltM a (A-59) 


NOTE: ASTROS does not support thermal gradient or moments s0 that the above 


equations are provided for completeness only. 
The in-plane stress vector (o)z at fiber distance z from the mid- 


surface is: 


oy Cx Cy 
lolz = Joy f = G- 2 poy f + G+ qoy (A-60) 
TxyJz Txy/1 TxyJ1 


where the stress vectors (oy ,eystxyl1? and (ox ,0ys%xy)27 are the bottom and 


top fiber stress vectors, respectively. 
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It a temperature Ty ts specified at the point where outer fiber 


Btresses are to be calculated, the additional thermal stress due to the 


difference between the specified temperature and a temperature that would be 
produced by a uniform thermal gradient T’ or thermal] moments (Mi_z is calculat- 
ed using: 

{body = [Go)lay}(Ty - Tp. - T's) (A-61) 
for a thermal gradient T’, and 


(M) 
{holy = cares ~ [Go}la)T, (A-52) 
A.7 FORCE RESULTANTS 
The forces at the wid-surface are evaluated by taking tre average 


stress values over the element thickness: 


(A) Forces 
Py 
(Fi = By ~ ({o}z1 + fo)z2) = (A-63) 
xy 
(8) Moments 
My 
(Yoo= 4 My = (lois, - (oy) 4 (A-64) 
Mxy 2 
(C) Transverse Shear Forces 
Qx 
(Qi = - Cirtey + (riz2) = (A-65) 
Qy 2 
¥ 


where stress vectors jojy] . (fo)y> are stresses at the integration points 


(default option) or at grid points (if requested) and, similarly, (r),, and 


ir),) are the transverse shear stresses. 
A.B THERnAL LOAD VECTOR 
The thermal ivad vector is computed as: 


(Pr) = fF (By [G}ied_ adv (A- 863 
Vv 


where the load vector (Py) Is defind «s: 


where (Fr) and (My) are the thermal forces and moments, respectively. 


The thermal strain vector is: 


«My oy 
(ey) = teeee = {2} AT (A-68) 
“Br ag 


where (*M)y and (€B)y are the thermal membrane and bending strains, and corre- 
spondingly {ay) and (ag) are the thermal coefficients of expansion for men- 
brane and bending. AT is dependent on the temperature loading being speci- 
fied. 


(A) For a specified grid point temperature the thermal membrane 


strain vector, (*M)r, is: 

fem) = fay)(Ty - To) (A-69) 
Ty = Grid point temperature 
To = Reference (stress-free) temperature 


(B) For a thermal gradient, the thermal bending strain vector, (€B)7, 
is: 


(ep) = (ag) oe T’) (A-70) 


(C) For thermal moments, the thermal bending strain vector, {€B)r, 
is: 


(ep} = [G2] (MI rn (A-71) 


NOTE: ASTROS does not support thermal gradients or moments so 


that the above equations are provided for completeness only. 
A.9 LAMINATED COMPOSITE MATERIALS 


The capability to model a stack of layers with a single QUAD4 element 
is detailed including the computation of equivalent "single layer" properties, 
i.e., membrane, bending transverse shear and membrane-bending coupling. The 
recovery of element forces, layer and interlaminar shear stresses and the 


computation of ply failure indices is also described in the following overview 
of theory. 
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A.9.3 


the 


OQvexview of Theory 


eslculation of the “"oversil" properties for the laminated ccmpos- 


ite elements fis based on the classical lamination theory with the following 


assumptions: 
(A) 
(B) 


{C) 


The 
in the folle 


WIE KE 


and the rema 


Each of the lamina is in # state of plane stress. 
The Laminate is preswaed to consist of perfectly bonded lamina. 


The bonds are presumed to be infinitesimally thin and non-shear 
deformable. That is, the displacements are continuous across the 
lamina boundaries so that lamina can not slip relative to one 
Thus, 


"special" properties. 


another. the laminate behaves as a egingle layer with 


material properties of laminated composite materials are reflected 


wing force-strain reletionship: 


2 T 
t Gy J t4G, | f) eu 7 eM 
es [--cc cee] Te---- eee 
- t’c, | 162 | O K - KT (A-72) 
L 0 | 0 f{ ts Gs 1 
Fy 
- Fy F Membrane Forces per unit length. 
Fy; 
My 
~ My 2 Bending moments per unit length. 
Bey 
Veo 
x 
- j Transverse shear forces per unit length. 
Vv 
y 


{ning terms have been defined previously. 
and Gy, terms are defined by the folloving: 


~ : f {Ge} dz 


S = x? (Gp) dz fAoTA) 
wm Af oe (Ge) ae 
: 


The limit on the integration are from the bottom surface to the top surface of 
the laminated composite. The elastielty matrix [Ge] hes the foilewing form 


for isotropic materials: 


pa Jee KE 0° 
lev? ley? 
[Cp = SYM =e ny) (A-74) 
ley 
G 
Cw ES (A-75) 
2(1+v) 


For orthotropic materisis, matrix [Gy] is: 


E E, 
eae Sern ay Ee 
l-vyv2q l-vyr2 


E 
{CE} = SYM —4— 0 | (A-76) 
l-vyy2 


G\2 


Equation A-73 may be rewritten as: 


G ~ Loc [6y,)% (2p - 2x. 
[644], qe ede PR ERD) 

ze Se ae 77) 
[G - L £ [G6 2x - 2p.1) A-77 
L ij), Sal [Gyy}h (2K + Spy ( 


NN : 2 
[Gry] ~ aA. oe fegyiR (a 2 


where [cy 4] is the reduced moduli matrix mveluated for each Jaina K after 


i 


transforming th lamina property matity from the fiber to the siement material 
£& P J 


aXE Ss . 

ke @ol 2p.qz ere the top and ooctowm Aisiances of lamina K from the 
Beoeetyic misdle plane of the Lamimace, as {liustraced in Figure A-G, and N is 
the mumber ef laminae (€or piles). Nere that the plies are mmbered werially 


etatting with Jo oat the bottom Laver. The bottom layer fs defiecd as the 


surtace with the largest -z walue in the @eleweni coordinate system TE the 


smear i ee Re a ME ET EET ee ee ae 


option to model membrane-only elements 4s axercised, matrices {Ga}, (G3), and 
{G,] are set to zero. 
If the user defined element axis in not coincident with the elemenc 


material axis, the user specified transformation angle @y, which references 


the element X-axis, is added to the layer orfentation ang.c. The property 


matrices [G1], [G2], and [G,] are chen transformed to the user defined element 


axis using the following equation: 


[Gg) = (BIT (eq) [V) (A-78) 
where 

[ con26 sin26 eostaindg 
{U}) «= | sin2s cos28 -cos@sinég (A-79) 

i -2sinfcos@ 2sinécosé cos%@-sin@¢ 
The transverse shear flexibility (G3] matrix is defined by: 

} S11 G12 
{A-09) 


[G3] = 
Gi2 $22 


and the corresponding matrix tcansformed into the user-defin’ 1 alenent  eordi 


nate ‘stem is given by: 


[CG] = (WJ? (Gy) [W] (A-81) 
where 
cos§ = sginé | 
{W) o» | (A-82) 
| -sin# cose | 


The derivation of the transverse whear Floxibiiity matcix [G3] for the 


laminate is considered xt. 


The mean value of; the transverse ahear wodulus, G, for the iamfiated 


composite is defined in terms of the transverse shear strain energy, U, 


th ough the depth as: 


a 2? 
RE si? eee tame” ee aly (A-83) 
2Gt 2 Giz} 


Sounique wean value of roarsverse shear gtrain is assumed to exlet for beth: 


we ca: and y-components of the element coordinate syscer, but for 


ease oof 


co cn AB ceraneecmer ion a a ae 


See 


a ee . Se . . 


dtecussion, only the evaluation of an wacoupled x-compoment of the shear 
moduli will be illustrated here. From zquation A-83, the mean value of 


transverse shear modulus is written fn the following form: 


an m E EF fp RK. oz (A-84) 
Cy, ye fel TUL (GD, 
B.S 


where G is an “overage’ transverse shear coefficient used by the element code 
and (Gy)q 2c the lecal shear coefficient for leyer 4. Yo evaluate Equation A- 
B4, it is necessary to obtain an expression for [vyy(z)]. This is accom- 
plisped by assuring that the x- and y-componente of stress are decoupled from 
one another. Unis assumption allows the desired equation to be deduced 


through an examination of a beam of unit cross-sectional width. 


a 


7 A 4 
a 
bene ie eh eee 


The equitibyism conditions im the hnxvizontal direction and for total 


MomMeak Ave; 


or So 
cowie md oe RY (A-85) 
bz cons 
aM 
Vy 4d aend om 0 (A-&6) 
ax 


Tf the location of tne meurrel surface is denoted by zy end p is the radius of 


curvatore of the beams, the ayvlel stress, o,, is expressed in the form: 


e x | 
iy os aa ora My. {A-87) 
(EE. 


AG 


_ SACL, OE COE TT TE EE SS SE SA ER TOA AEE at SER FAA RM CS EEK. LAR: 
Aes P mn 


Equation A-87 is differentleted with rempect te x and coubined with Equations 
A-BS and A-86. Fox constant E,, the result is intsgrated to yleld the follow. 
ing expression: 

v Pr 


Zz 
om One eed. dz ew ] Ex We CLC 
i (Ey x 2 | X 4-3 L 


Tee (A> 88) 
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Equation A-88 fs used in the esalysis of n-piy Laminates because 


sufficient conditions exist to determine che cometante (, Ciet,2,.....n) and 


the “divectional bending center,” zg. Por exemple, wonsider cho Following 


laminated configuration: 


RW 5-2 RENEE OE SIE SOS 0 I aca | 


* WIPO ELS EETIOR. 


Bors mance anes sone nfm SS 


a ® fan ft 
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At the bottom surface (1-1, e-z,, and ry,70), therefore: 


ee 
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a zee 
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Slack ply at the interface between plies fel and fe? Crezy)s 
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Then, in the ply, 2] < Z < 22, the shear is: 
VLE - 9 2 
e rx2(zZ) = (7xz)y + Bek? [zy (z-21) - 2 (22-21) (A-92) 
(EI)y 2 
In general, for any ply z1., < z < zy, the shear is: 
VLE - 2 2 
bd teo(Z)g = (1 xz) gq.) + RAE [2 (2-24.21) - 2 (2?-24.4)) (A-93) 
(El), 2 
At any ply interface, z;, the shear is therefore: 
Vv i - 3 
@ Tx - —t £ T - + 24.3) (A-94) 
CRS Ge gn eg a 
where Ts = Z4-Z4-1- 
Note that the shear at the top face, (%xz),, is zero and therefore: 
@ 
Vv fj. on n (z,4z ) 
(xz), = —%— fz, Ty = T, —J-i:1'| - 0 (A-95) 
7 m (EI), le j=l Exj J jel Exj J 2 


Equation A-95 proves that if Z, is the bending center, the shear at the top 


surface must be zero. 


A better form of Equation A-93, for this purpose, is: 


VLE - 2 
e [?xz(2)}, = an fxy + z(z~-z4.1) - : (22-244) | (A-96) 
x . 
where 
1-1 a 
© fxi = Be, re Fxj Ty [2x - 2 (2ytzy-1)] (A-97) 


Substituting Equation A-96 into Equation A-84 yields: 


n 
1.» IL. y 1k, , (A-98a) 
bd Gy (E1)2 J=1 Cx, | 

x 


5] 


where 


Rep = (Exi)? ter + (2y°24-)Ty - 4 15?) fx4 + Azy-224-1) 


a3 3 3.3 
- r Ty) zyTq + a 21-4 + 7 24.17, + rs Ty)Ty | (A-98b) 


This expression for the inverse shear modulus for the x-direction is general- 
ized to provide for the calculation of each term in the two-by-two matrix of 


shear moduli as: 


as = n g ol “i 
[Cyg]) = ——_ 5 [C,}) Rui (A-99) 
(EI)2 isl 
kk 
where 
k = 1,2 
1 - 1,2 


Note that if no shear is given, tolj-t.0, and also that, in Equation A-99: 
= * 
(EI)}] = 1,1 term of I x [G2] 


— * 
(EI)92 = 2,2 term of I x [G9] 


where [G*9] is calculated in the game manner as [G3} except that Poisson's 
ratio is set to zero. The moduli for individual plies are provided through 
user input. Because Gj)9*G21, in general, an average value is used for the 


coupling terms. 


c 


G11 12) ave 
[G3] - le (A-100) 
(712) ave G22 
A.9.2 Element Layer Stress Recovery 
The linear strain variation is given by: 
(ex) = (ey) - 2(K) (A-101) 


uy 
np 


where 


“ayer strain vector in the element coordinate syatem. 
Reference surface strain in the element coordinate system. 


Reference surface curvatures in the alement coordinate 
aystem. 


Distance of the mid-surface of the layer k from the laminate 
reference surface. 


The individual layer stress vector in the fiber coordinate system Is: 


where 


{oy} 


oy) 


{Gy] 


[T) 


(ey) 


(Gy) {T] 4#y) (4-102) 


Layer stress vector in the fiber coordinate systen. 
Stress-strain matrix in the fiber coordinate systen. 


Transformation matrix to transferm strains from element 
coordinate system to fiber coordinate system. 


Layer strain vector in the element coordinate system. 


For element teuperature and/or thermal gradients, the strain vector 


has to be corrected for thermal effects before applying Equation A-103: 


(ey 


) 


ley} ~ (a) (TP + 2T') (A-103) 


and for thermal moments 


where 


lex 


, 
léy 


{a} 


) 


) 


Cay 


(ex) = (e)7 (A-104) 


Mechanical gtrains. 


Thermal coefficients of expansion in the elemert coordinate 
systen. 


Element temperature. 
Element thermal gradiert. 


Listance from the middle of the layer to the jJaminate 
refcorence surface. 


Leyoy stvains due to thermal momenis ip the element coord) - 


nate system. 


eee 


The thermal strain vector due to applied thermal moments is determined 
| by substituting for (M} in Equation A-73 and solving for the reference surface 


strains and curvatures, (eT) and {KT}, respectively. 


A.9.3 Interlaminar Shear Stresses 


The interlaminar shear stress Tyz, Txz can be computed at any ply 


interface from Equation A-96. 
A.9.4 Force Resultants 
Forces and moments for the element are computed using: 
N 


(F) = JE, lox) Ty 


i-1l, N (No. of layers) (A-105) 


N 
{M} = Fa -24Ty (oy) 


where 
{F} = In-plane force resultants. 
(MX) = Qut-of-plane moments. 


Stresses in the element coordinate system. 


~ 
Q 
bad 
! 


Ty ~- Layer thickness. 


Zz ~ Distance from the middle of the layer to the laminate 
reference surface. 


A.9.5  Fallure Indices 


Failure indices assume a value of one on the periphery of a failure 
surface in stress space. If the failure index is less than one, the lamina 
stress is interior to the periphery of the failure surface and the lamina is 
assumed “safe” and if it is greater than one the lamina is assumed to have 
"failed." The failure indices represent a phenomenological failure criterion, 


because only the occurrence of failure is predicted. 


The analytical definition of a failure surface in stress space for a 
lamina subjected to biaxial (planar) states of stress is provided via the 
following failure theories. 
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(1) HILL 
(2) HOFFMAN 
(3) TSAI-WU 
(4) MAXIMUM STRESS 
(5) MAXIMUM STRAIN 


In the analysis of laminated composites, which are typically orthotropic 
materials (possibly exhibiting unequal properties in tension and compressicn), 
the strength of orthotropic lamina is a function of body orientation relative 
to the imposed stress. In order to determine the structural integrity oY the 
launina, a set of intrinsic strength properties (allowable stresses or allowa- 


ble strains) in the principal material directions are defined as: 
Veo. Ultimate unlaxial tensile strength in the fiber direction, 


Od Wlcimate uniaxial compressive strength in the fiber direc- 
tion, 


Ye = Ultimate uniaxial tensile strength perpendicular to the 
fiber direction, 


Yo = Ultimate uniaxial compressive strength perpendicular to the 
fiber direction, 


S - Vitimate rlanar shear atrength under pure shear luading, 
Ey ws Ultimate uniaxial tensile strain in the fiber direction, 
Ee Ultimate uniaxial compressive rtrain in the fiber direction, 


Feo Ultimate untaxial tensile strain perpendicular to the fiber 
direction, 


Feo Ultimate: uniaxial compressive strain perpendicuisr to the 
fiber direction, and 


E; oo Ultimate planar shear strain under pure shear loading. 


Fur most composite materials, the planer shear strengths and strains 


are equal for positive and negative whee: loadings. 


The five fallure ctheerfes und m bonding Eallure Index are now de- 


ecribed: 


ff, 


HALLS THEORY 
of of Cod r2 
wh + he - ah 24 12 & FAILURE INDEX (FI) €A-106) 
x22 x2 $2 


anc Zeek, if o, fis positive, and KX, if o) iw negative; eflailarly for y. For 
the interaction term, (9192) 7x2, KeX; Lf oy79 in positive K=X, etherwisze. 


HOFFHAN'S THEORY 


Note that this theory takes into account the difference in the tensile 


and compressive allowable stresses by usirg linear terms in the failure 


equation, 


TSAL-WU_THEORY 
This quadratic interaction theory allows for the strength predictions 
wherein interaction mnong stress components can be considered in determining 


strengths in a biaxial ffeld. Thus, in the case of an orthotropic lamina in a 
general state of planar «tress: 


9 


2 2 
Fyoy + Foo2 + Fy,y0) + Fonog + 2Fyi 20302 + Fegrig = FI (A-10€) 
Fp -e bh e-k Fo ~ 2-1, 
XX, Ye Ye 
fA-109) 
Fi ae Fone ae » %-F? ae pane ; Fg6 ~ pe 
XrX5 ¥eYe $ 
However, 


and F)9 needs to be determined experimentally, from » biaxte’ test. 


satisfectory resulce pay be obtained by setving it te zero. 


Fallure 48 sasumed toa occur when any one of th «tvess components fs 


equal to {ts corréspendiag Gstrinsic strangth property. In mathematical roru, 


the Maximum Stress theory ts given br: 


oy 2 ky, of > 9 + OF S -Ne, ay < 0 
one Ye, 02 > 9 3 of S Ys, op <0 (A-110) 
r92 2S, eyy > 0 5 wig SS, #42 <G 

where the intrinsic strength propercies are as defined previously. 


MOXIMUM STRATK 


dhe Maximum Strain theory fe enalegous co the Maxiccem Streas theory. 
Faflure is assumed to result when any one of the strain componenta is squal to 
its corresponding intrinsic ultimate atrain. In mathematica’ form the Maximun 
Strain theory is given by: 


¢y # Ec, «7 > 0 | « 
eg 2 Fe, 29 >0 3; eg S Fo, €2 < 0 “A-111) 
12 2 Eg, vi2 > OC F Wye 8 By, ya2 < 0 
where the Intrinsic ultimate etrains are as defined previously. 
FSILVRE TNEEX OF PONDING 
The Kc<alluce dncex of bouding material ds calculated ts che maxon 
interlaminay shear etreas divided by the allowable bonuing rtress. 
A.10 COMRECTION OF OUT-OF- PLANE SHEAR STRAIN 


The typical foraulation for # QUAD4 type finite elament fellows a 
standard bilineay iaoparnmetric theory, wath dlvertionsl reduced integration 
for out-of-plane shear strain, However, this formulation ues been faund se be 
inadequate when the geometry ot the element fs drregniear, and # correction 


vefined heceln bas been implemented In ASTROS te correst this preblen, 


The wodification is based wpon the theery preseutec by Hughes and 
Tezdayax (Referei te A-1), but ds generalized tu inciude non-pianarity of the 
element, and special featurer te accommodate ASTROS: s atructure. The flormula- 
tion enforces constant swbear along wach edpe of the erecent, alinminating che 


need to perform reduced integretfcn. 


The form-Tation of this medification co.u.lets of establishing otras. 
{> 


dlaplacement xelacionshivs in the clament coordinate system. Te invo 
degrees of freedom (dof). the retetienal part of which will be modiffed later 


to include the singular.ty about the normal to the mid- su. face 


A.10.1 Geometric Varisbles 


The following terms are defined for each 


edge of an irregular-shaped, 
non-planar element: 


A Unit Normal Vector ( © ) in the direc- 


tion of the next node as {llustrated; 
6 


1 @ 2 
. 1 
A Unit Normal Vector ( T ) which is a 
normalized average of the nodal normals 


to the mid-surface along that edge; 


Length of each edge (hy) 


; and cosine of the internal angle at each 
corner (a;). 


A.10.2 e 


Given the following numbering sequence. 


At the middle of each edge, the constant 


shears parallel to edges a, b, c and 
d, respectively, are: 
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gE 7 oe Ty, > (Up - Ug) - de, (8p + #a) 


ae - (Ua - es . rv} ? 
B,. he m, * Wa - Up) 2 eb fla + 9p) 
(A-112) 


o - ~ ~~ 
or A me + (Ug - Bed - & (Oa  8e) 


~ ~*~ =» “y ~s ~ 
B, * dmg + We + Ug) - dog Ue + a? 


es + 
where uw ana @ are the vectors of translations and rotations at each node, 


respectively, in the element courdinete syster. 


The shear vector (7) at node (6) is given by: 


mh = — (gp. + Bath)? # duck Be + 2 ai en {A-113) 
ee " eae 6 
or 
Yb “ (eat ey) (iy ay. aI (egrapey) (Da Up) 
Q-a? 2 ha a2 dy 


- ae Or rere ay) . chat ce timp i 


fee ‘3? yhp Qa of yh 
(A-1ih) 
aoe a ate) 
- [ann (e, age dest $9} 
2(4-0°5 ; 
& ; 
r 7 f+ ast ad +e 
oe aes, arte, : inate, awuels en: aa 
beh oe as 


= = 


ae ey ton ecient Bi, ) 


Cowie semen 


ar 


and similarly for the other -des, by permtations of the a, b, c¢ and a 


subscripts. 


A.10.3 Nodal Contributions of Shear Strain 


The contribution of each node to the total shear strain (77) evaluated 
at an integration point is: 
- 4 ~ 
wo 2, Nav (A-115) 
: deol 
The "pseudo-contribution" of each edge to the tetal shear strain (G) has the 
follewing form: 


~ N ~ ~~ N red 
GC. = -—-& (eg + ageo) + bh. (ea + apep) 
: a 1-a2 1-a’ 
a a 
3 N ~ ~ N aa r 
G = RB. (ey + apes) + —S. (oy + ageg) 
b leat 1-a2 
b 
(A-116) 
a N + os W > ~~ 
Go —- fe, + apeg) + —H. (e, + ageg) 
© L-a? l-a? 
c & 
é Boe (entered ae emece.) 
: Wace ae €., + Aro le Sa fd + Bele 
d l-aé lea? 
a c 


Hence, the colwmns of the |B) matrix partition for shear, corresponding to 


node b, [BS,], are: 


an ni - ni> t 
{Bony} = —B Gy - ob Ch | 
ha Ob i {1.2.3 
r (A-117) 
at _ get ne j fet 58 


A.10.4 Jronsformaiious 
The follewing transforaactons heve to be porcormed before the preced 
4ng formulation ean replace the @asiseving |—§5) mwateds wenerarfion for out-of, 


plane shear. 


where 


> ¢ 


= 14) ) @ 
(TIE), (BSy].) [ooo ee . (A-118) 


Is the orthogonal transformation between integration 
points and the element coordinate tys em, required 
since a1] the strains ere calculated in the 1 aystiu. 
Is a 3x3 identity matrix. 


Is the 3x3 transformation which taxes into accourt the 
following facts 


«A; Hughes’ convention for rotations is different then 
the one implemented in STROS; and, 


(B) The rotation about the werual to the wid-surfare 
at each grid point is elagular. 


Tf NV is the normal veers ut a given grid point, chen. 


Ci wee wre y 


[TEE]: = nv3 9 -pyh (A-119) 


Figure A-1. Internal Element Coordinate Syscem 
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Figure A-3. Deformations at Grid Point i 


Figure A-G. Deformations in the Gisbal Otrection 
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Tigure A-S5. Material and User Defined Element Axes 
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Figure A-6, Geometry of @ R-“aye ed Element 
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APORNDIX a: GUEDSOUTNES POR MODELING WiTa OU: 
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GUIDELINES FOR MODELING WITH QUAD4 (TRIAS; ELEMENYS 


ASTROS (an Automated STRuctural Optimization Syatem) and four versions 
of NASTRAN (COSMIC NASTRAN, UAI NASTRAN, CSA NASTRAN and MSC NASTRAN) provicie 
QUAD4 and TRIA3 elements for the analysis of plates. They are basically fiat 
plate elements, but they are of.en used for shell etructures as well as 
approximations. These approximations become closer to reality with a finer 
mesh size. Much of the discussion in these quidelines is with reference ts 
the QUAD4 element, but comments are equally applicable to the TRIA3 element. 
the original QUAD4 formulation in programs ASTROS, COSMIC NASTRAN aid UAT 
NASTRAN ig just about identical, and ditferences (if any) ere the result of 
subsequent revisions. MSC NASTRAN and CSA NASTHAN are different, but never- 
theless all five programs give comparable results for most well posed 
problems. Also the input requirements and the card structure are guite 
compatibie in all five cases with minor differences in interpretetion. The 
purpose of these seminar notes is to explain various QUAD4 modeling options 
and to point-out the differences in the five progres. An explanation of the 
input parameters on the various cardi: is co:.sidered che easinss way co 
accomplish this objective. 


The QUAD4 is one of the most extensively used elements in Ne! “RAN as 
well as ASTROS. It is a very versatile element and can be used to model a 
variety of plate simulations such as (see Fic. 1) 


Membrane (inplane lcading} behavior 

Bending (out of plane loading) behavior 
Membrane-bending (uncoupled) 

Membrane-bending (coupled-linear) 

Laminated plates 

Layered composites 

Sandwich plates with metal face she» s 

Sandwich plates with layered composite fave si. ets 
Isotropic materials 

Anisotropic (including orthotropic) materials 


Ur Tome aners 


Application of the QUAL4 is often confusing becaus of the many optics 
available for its use. 


There are five cards which describe the inp. yuraneters for the QUADS. 
They describe its geometry and properties along wath some avxiliary informa- 
tion, 


Geometry and Property Cards 


CQUAD4 - Connection card 
PSHELIL - Property card for homogeneous and randwich platen 
PCOMP 
or 
PCOMP1* - Property cards for Jaminace. layered plates 
or 


PCOMP2 * 


Material Cards 


MATL ~ Isotropic materials 

MAT2 - Anisotronic materials 

MATB - Orthotropic materials 

FLOAND4 - Pressure load definition on the QUAD4 element 


The property cards PCOMP1 and PCOMP2 are not applicable to MSC NASTRAN. 
For a given element either the PSHELL or PCOMP card is applicable but not 
both. PSHELL cards are for homogeneous ({nonlaminated) and sandwich plates 
with nonlayered face sheets. PCOMP cards are for laminated (layered) plates. 
In the case of sandwich plates with layered face sheete the honeycomb 
(sandwich) core will be treated as a laminate or layer. 


A supplementiry explanation of the parameters on each of these cards 
should aid in understanding the modeling nuances of the element. 


Notes on the CQUAD4 


The format of the CQUAD4 is the same for all the NASTRANS and ASTROS 
with the exception of an additional parameter "TMAX" (field 3 on the 
continuation caid) in ASTROS. "TMAX" is the maximum allowable thickness of 
the plate, applicable only in optimization. 


The definitions of the parameters in fields 2 to 7 are self explanatory 
and need no further clarification. Similarly no additional explanation is 
necessary for the thickness parameters specified in fields 4 to 7 on the 
continuation card. However, the parameters TM (THETA) and ZOFF need a 
supplementary explanation or caution. 


Parameter TM (THETA) 


The parameter TM defines the material property orientation. There are 
two options for this definition. 


Option 1: 


Define the angie between the side of the element (connecting G1 and G2) 
and the material axis. This is the least desirable opticn. It is vone to 
errors, because every time the sequence of the element ccnnection changes, 
the angle must be changed. Also in a complex three dimensional motel it is 
not easy to determine this angle without writing a prey rocessor. 


Option 2: 

The integer option is preferable. An integer i:: field 8 refers to a 
separate coordinate system for defining the orientation of the material axis 
of the element. The materia] property definition is now independent of the 


connection sequence. The new coordinate system : an be defined with a CORDZR 
card. 


The offset parameter provision in the QUADA elemer: constitutes a siani- 
ficant enhancement fer plate elements. Refere the GUsba the grid points of 
the structure conuid only be defines on the mid-surtace of the plate elements. 
The Bar (beat ot bend! was the only other element wit’ an offser capability. 


Heowever 


, &ome of the mass elements have the offset capability. 


The offset, ZO, is shown for various cases in Fig. 3. Note the 
aistinction between the grid point surface and mid-surface of the element. 
The definition of the offset ZO on the CQUAD4 is the same in all four 
NASTRANS and ASTROS. The offset implementation in COSMIC NASTRAN and ASTROS 
appears to give more consistent results than MSC and CSA NASTRAN. See 
Reference Al for results on the benchmarking of offsets. 


Notes _ on PSHELL 


PSHELL or PCOMP are the property cards referenced on the CQUAD4 (in 
field 3). The PSHELL is to be used when the plate is not laminated (or 
layered), while the PCOMP is for laminated plates. Only one of these is 
applicable for a given element. A preprocessor in NASTRAN (ASTROS) generates 
equivalent PSHELL cards from the PCOMP cards before proceeding to the solu- 
tion. NASTRAN (ASTROS) provides versatility to the QUAD4 element through the 
PSHELL card. It allows the modeling of membrane (inplane), bending, shear 
and membrane-bending coupling behavior. Fig. 2 illustrates key features of 
the elements described on the PSHELL card. It is a sandwich plate with two 
face sheets separated by a honeycomb core. 


The first two fields of the PSHELL card are for the name and property 
identification called from the CQUAD4. The third field, MID1, is the 
material identification number for the face sheets in membrane behavior. The 
parameter T is the total thickness of the two face sheets. MID2 is the 
material identification number for bending behavior, MID3 for shear and MID4 
for membrane-bending coupling. There are two types of membrane bending 
coupling. The coupling resulting from asymmetry in plate construction (non- 
symmetric laminates) is called linear coupling. Nonlinear coupling, on the 
other hand, is a result of the interaction of internal forces such as inplane 
and out of plane (beam-column effect) forces. The latter coupling can be 
accounted for only in differential stiffness and/or buckling analysis. The 
parameter 121/T? (field 6) can be calculated by using the following defini- 
tion for I. 


I is basically the moment of inertia of the face sheets about the neutral 
axis (centroidal). It is assumed that the face sheets are symmetric about 
the neutral axis. If they are not, the moment of inertia about the neutral 
axis can be calculated. For solid plates this parameter is simply 1.0 


The definition of the parameter TS/T is obvious from the figure. 


No further explanation is necessary for the parameters in the next three 
fielcés, NSM, Z1 and 22. 


The parameters MCSID and SCSID (not available in MSC and CSA NASTRAN) 
refer to the material coordinate system. As stated in the description of the 
card CQUAD4 there are two options for this definition. By leaving the field 
blank or a real value the first option is invoked. In this option the 
parameter represents the angle between the side of the element connecting the 
grid points Gl and G2 and the material axis. The second option is an integer 
which refers to a coordinate system defined on a COORD-card. The second 
option is the most desirable because the grid point sequence on the CQUAD4 
card does not affect the material axis. 
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The offset parameter ZO is the same as defined on the CQUAD4 (See Fig. 
3). An important point to note is that there is no provision for the offset 
definition on the PSHELL card in MSC and CSA NASTRAN. There is some 
advantage in having this option on the PSHELL, because when the number of 
PSHELL cards is significantly fewer than the number of CQUAD4 cards, this 
parameter need not be repeated on all the CQUAD4 cards. The entry on CQUAD4, 
however, overrides that on the PSHELL card. 


The PSHELL card provides the facility to model homogeneous as well as 
sandwich plates. However, the face sheets of the sandwich plates are assumed 
to be homogeneous (isotropic, orthotropic or anisotropic) plates. A 
discussion of sandwich plates with face sheets made of layered composites is 
differed until the PCOMP cards description. The minimum thickness parameter 
TMIN is applicable only in ASTROS optimization. 


Before leaving this discussion, it is worth pointing out some anomalies 
(idiosyncrasies) (as they exist at present, Feb 1993) in the application of 
the PSHELL card in various versions of NASTRAN (ASTROS). 


Specifying MID1 and leaving the remaining three material identifications 
blank simulates the membrane behavior in all five programs with little or no 
differences to point out. 


Specifying MID2 only invokes bending behavior only (as it should be) in 
MSC NASTRAN. It does not compute shear deformation. CSA NASTRAN on the 
other hand includes shear deformation, ostensibly with the material proper- 
ties invoked by MID2. ASTROS and COSMIC NASTRAN go even further by including 
membrane, bending and shear deformations even though only MID2 is specified. 
Membrane and bending behavior are computed by the material properties called 
from MID2. This is not so bad, -because these two behaviors are generally 
uncoupled. The most important point to note is that ASTROS (COSMIC NASTRAN) 
computes the shear deformation by assuming a material infinitely stiff in 
transverse shear when the MID3 field is blank. These results are, in 
general, not acceptable. The easiest way to avoid shear stiffness over- 
estimation at present is not to leave MID3 blank when MID2 is specified. See 
results of the examples at the end of this Appendix. Future versions of 
ASTROS and COSMIC NASTRAN will uncouple these computations. 


Notes on PCOMP, PCOMP1, PCOMP2 


The purpose of PCOMP, PCOMP1 and PCOMP2 is to define element property 
parameters in modeling laminated plates including layered (fiber reinforced) 
composites. All three cards serve the same purpose except the options are 
different. If the layers are made of different materials and the thicknesses 
of the layers are all different, then the PCOMP card is appropriate. If all 
the layers are made of the same material and thickness, then PCOMP1 is 
appropriate. If the material is the same, but the thicknesses are different, 
then PCOMP2 is appropriate. The first two fields on the PCOMP cards need no 
further explanation. 


Parameter 2@ 


The parameter Z@ refers to the distance from the grid point surface to 
the bottom of the plate. The plate bottom surface is defined in Fig. 4. It 
is the reference surface from which the stacking sequence of the laminates 
is defined. The parameters ZOFF defined on the (QUAD4 and PSHELL) are a 
source of confusion sometimes. However, figures 3 and 4 should provide 
better clarification. 
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Parameter SBOND 


The bonding material shear stress is indirectly related to the 
interlaminar shear and its value is generally empirical. A value of 400 to 
500 psi for SBOND appears to be reasonable in the absence of a value obtained 
from experiments. Any approximation of this parameter will not affect the 
analysis results. It affects only the failure theory which is basically a 
postprocessing function. 


Parameters FT (Field 6) 


This parameter simply identifies the desired failure theory. The five 
failure theories are discussed in detai] in the theoretical secticn of the 
QUAD4 seminar notes. The possible values ot FT are: 


"HILL" for the Hill Theory 

"HOFF" for the Hoffman Theory 

"TSAI" for the Tsai-Wu Theory 

"STRESS" for the Maximum Stress Theory 
"STRAIN" for the Maximum Strain Theory 


Parameter LOPT 


This parameter description needs no further explanation than what is 
given in the PCOMP card descriptions. It should be pointed out, however, 
that no parameters are in fields 7 and 8 in COSMIC NASTRAN, while MSC and CSA 
NASTRAN define "TREF," a reference temperature for thermal stress analysis 
and "GE" the material damping coefficient for dynamic analysis. 


Parameters MIDi, Ti, THi, and SOUTi 


These parameters pertain to the ith layer. MIDI is the material 
identification number of the first layer. The layer count goes up from the 
bottom surface of the plate. For the definition of the bottom surface see 
Fig. 4 for the definition of 26. The MIDi refers to one of three material 
cards: MAT1 for isotropic materials, MAT2 for anisotropic materials and MATS 
for orthotropic materials. The parameter Tl defines the thickness of the 
first layer and TH1 refers to the orientation of the material axis with 
reference to the material axis defined on the CQUAD4. SOUT1 is the stress 
output parameter. Then the parameters are repeated for all the layers unless 
the symmetry option is used under the parameter LOPT. If any MIDi, Ti or THi 
are blank, then the last non-blank values specified for each will be used. 


Material Cards 


Isotropic & Anisotropic 


The parameters on MAT1 and MAT2 are self explanatory from the card 
descriptions. 


Orthotropic Material MATS8 


Most of the parameters on MATS are self explanatory with the exceptions 
of G1Z and G2zZ (fields 7 and 8). When these parameters are left blank, MSC 
NASTRAN and CSA NASTRAN do not calculate the transverse shear deformation in 
plate bending problems. ASTROS and COSMIC NASTRAN, on the other hand, assume 
that the material is infinitely stiff in transverse shear, and thus over- 
estimate the stiffness of the element. To avoid such overestimation and to 
obtain comparable results to the other NASTRANS, transverse shear values have 
to be provided. Values of about two or three orders of magnitude less than 
the modulus of elasticity of the material are recommended. 
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These guidelines and the accomparying examrrlss should clarify most of 


che questions arising in the use of the QUADa and TRIA3 elenenta. 
REFERENCE 


1. Pitrof, Stephen M. anc Venkayya, Vinoperlsa B., Benchmarking the 
QUAD4/TRIA3 Element, Twenty-first N&sTuAN liser’s Colloquium, Tampa FL, 
April 1993. 
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INPUT IS PROVIDED 


PROBLEM #5 


RECTANG!AR PLATE 


FIGURE 5 


A finite element model of 1/4 of a rectangular plate is shown in Fig. 5. The 
length of the plate is 2 or 10, the width is 2, i.e. Aspect Ratios (AR) 1.0 
and 5.0, and the thickness is .001. The material properties are given as 
E=1.7472x10 and v0.3. The plate is subjected to two loading conditions 

and boundary conditions for each aspect ratio as follows: 


CASE 1: Clamped Supports. Concentrated load of P = -4.0x1074 at the center 


of the plate. 
CASE 2: Simple Supports. Uniform pressure load of q=1074 over the plate. 
Convergence is studied by varying the mesh size. Input data is given for a 


6x6 mesh. Theoretical results for the lateral displacement at the center of 
the plate are given in Table 2. 


AR = 1.0 
223 


TABLE 2 
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Example A4: The plate is modeled as a sandwich plate with composite face 
sheets and a honeycomb core. The thickness or the core is 2in. 
Two cases are considered. In Case 1 the MATS entries Glz and 
G2Z are not specified. In Case 2 G1lZ and G2Z are given the 


value 4.0X10°psi. 
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EXAMPLE A3: QUAD4 - SHEAR - G1Z, G2Z SPECIFIED 
Composite Solid Plate 


Z DISPLACEMENT (x10) AT THE CENTER OF THE PLATE 


EXAMPLE A4: QUAD4 - SANDWICH - COMPOSITE FACE SHEETS 
Composite Sandwich Plate with Honeycomb Core 
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APPENDIX B: SAMPLE PROBLEMS 
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PROBLEM 
NUMBER 


1° 
2 


3 


12 


12 
13 


14 


QUAD4 SEMINAR 


APPENDIX B 


SAMPLE PROBLEMS 


DESCRIPTION 

PATCH TEST FOR PLATES 

STRAIGHT CANTILEVER BEAM 

CURVED BEAM 

TWISTED BEAM 

RECTANGULAR PLATE 

SCORDELIS-LO ROOF 

SPHERICAL SHELL 

COMPOSITE PLATE - PURE TWIST LOADING 
COMPOSITE PLATE - UNIFORM PRESSURE LOADING 
OPEN COMPOSITE TUBE 

STRAIGHT TEAM TEST - SIMULATION OF EQUIVALENT 
ISOTROPIC PROPERTIES. LAMINATE CONFIGURATION 
0/0/0/0. 

COMPOSITE SHELL ROOF MODEL 

COMPOSITE RECTANGULAR PLATE 


COMPOSITE SANDWICH PLATE 


PAGE 
NUMBER 


83 
85 
86 
87 
8B 
89 
90. 
91 
92 


93 


94 
96 
97 


99 


* Problems 1 through 7 reference isotropic material and are defined in the 
MSC/NASTRAN APPLICATION MANUAL, SECTION 5. 


+ Problems 8 through 12 are provided by UAI/NASTRAN. 
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B. 


1 


PROBLEM NESCRIPTIONS 


INPUT IS PROVIDED 


PROBLEM #1 


PATCH TEST 


FIGURE 1 


A model of a patch test for plates is shown in Fig. 1. The length of the 
plate is .24, the width is .12, and the thickness is .001. The location of 
the inner nodes is given in Table 1. The material properties are given as 
E=1.0xl0° and v=0.25. Displacement boundary conditions are specified for 
the following cases: 


CASE 1. Membrane Plate 


" 


u = 1079(x+y/2) 


1073(y+x/2) 


it} 


v 
CASE 2. Bending Plate 


W = 1073(x2+xyry?) /2 


@ = W = 4979 
3y 10° (y+x/2) 


x 
9 = . OW = 197 3(-x- 
y 3x 10° ~(-x-y/2) 
The theoretical solution for Case 1 is given by 
Ey ey Y 10 ~;3 Oe oy = 1333.; Ty = 400. 


The theoretical solution for Case 2 is given by 
Bending moments per unit length: 
- 2 =7 z 
m= Mm = ].111x10 My = 10 
Surface stresses: 


Roepe +.667; Ty * +.200 


TABLE 1 
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INPUT 1S PROVIDED 
PROBLEM #2 


STRAIGHT CANTILEVERED BEAM 


6 8 
a) Reguiar Shape Clements 


b) Trapezoidal Shape Elements 


c) Paratlelogras Shape Elements 


FIGURE 2 


A finite element model of a straight cantilevered beam with rectangular, 
trapezoidal and parallelogram shape elements is shown in Fig. 2 (plane view). 
The length of the beam is 6, the width is 0.2 and the depth is 0.1. The beam 
has material properties E=1.0x10° and v=.30. All the elements have equal 
volume, The beam is subjected to 4 loading conditiors as follows: 


CASE 1: A unit force is applied at the tip in the X direction. (Extension) 


CASE 2: A unit force is applied at the tip in the Y direction. 
(In-plane shear) 


CASE 3: A unit force is applied at the tip in the Z direction. 
(Out-of-plane shear) 


CASE 4: A unit force is applied at the tip in the +Z direction. (Twist) 


The theoretical solution for the displacement at the tip in the direction of 
the load for Cases is 3.0x10°~, for Case 2 0.1081, for Case 3 0.4321 and for 
Case 4 .03208x10 “. 
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INPUT IS PROVIDED 


PROBLEM #3 
CURVED BEAM 


90° 


FIXED 


FIGURE 3 


A finite element model of a curved cantilevered beam is shown in fig. 3s 
The inner radius is 4.12, the outer radius is 4.32, the arc js 90°, and the 
thickness is 0.1. The beam has material properties E=1.0xl10° and v=0.25. 
The beam is Subjected to two loading conditions as follows: 


CASE 1: A unit force is applied at the tip in the +Y direction, 
(In-plane shear) 


CASE 2: A unit force is applied at the tip in the +Z direction. 
(Out-of-plane shear) 


The theoretical solution for the deflection at the tip in the direction of the 
load for Case 1 is .08734 and for Case 2 is .5022. 


86 


INPUT 2S PROVIDED 
PROBLEM #4 


TWISTED BEAM 


FIGURE 4 


A finite element model of a twisted beam is shown in Fig. 4. The lenath of 
the beam is 12, the width is 1.1 and the depth is .32. The beam is twisted 
50° from root to tip. The material properties are given as E=29.0x10° and 
v=0).22, and the beam is subjected to two loading conditions as follows: 


CASE 1: A unit force is applied at the tip in the +X direction. 
(In-plane shear) 


CASE 2: A unit force is applied at the tip in the -Y dtrection, 
(Out-of-plane shear) 


The theoretical solution for the deflection at the tip in the direction of the 
load for Case 1 is .005424 and for Case 2 15 .001754. 


INPUT TS PROVIDED 


PROBLEM #5 
RECTANGULAR PLATE 


FIGURE 5 


A finite element mode! of 1/4 of a rectangular plate is shown in Fig. 5. The 
length of the plate is 2 or 10, the width +s 2, i.e. Aspect Ratios (AR) 1.0 
and 5.0, and the thickness is .001. The material properties are given as 
E=1.74?2x10 and wO0.3. The plate is subjected to two loading canditions 

and boundary conditions for each aspect ratio as follows: 


CASE 1: Clamped Supports. Concentrated load of P = -4.0x1074 at the center 
ot the plate. 


CASE 2: Simple Supports. Uniform pressure load of Gio” over the plate. 


Convergence 18 studied vy varying the mesh Size. Input data 1s given for a 
6x6 mesh. Theoretical results for the lateral displacement af the center of 


fal 


the plate are given an Table ?,. 


INPUT IS PROVIDED 


PROBLEM 6 
SCORDELIS - LO ROOF 


FIGURE 6 


A model of a Scordelis-Lo Roof i.e. a singly curved shell, is shown in Fig. 6. 
The roof has a radius of 25', a length of qo!» end.a thickness of 0.25'. The 
material properties are given as E=4.32x10° 1bs/ft™ and v=0.0. The roof is 
loaded by its own weight, and the weight of the roof is 906 Tbs/ft". Only 1/4 
of the roof is modeled. Symmetry boundary conditions are imposed on the 
interior edges and U,=U_=9 on the curved edges. Convergence is studied by 
varying the mesh sizé. *tnput data is given for an &x? mesh, The theoretical 
sclution for the midside vertical displacement is .3086'. A value cf .3024' 
was used for normalization of the results. 


Reference: Zienkiewicz, O.C., The Finite Element Method, Third Edition, 


McGraw-Hitl Book Co, (uk) Limited, pp 349-351, 


P 


INPUT AND OUTPUT ARE PROVIDED 


ROBLEM #7 
SPHERICAL SHELL 


yvF @ 2,9 
(on quadrant) 


F © 2.0 (on quadrant) 


FIGURE 7 


A model of a spherical shell, i.e. a doubly-curved shell, is shewn in Fig. 7. 
The shel’ has a radius gf 10 and a thickness of 0.4. The material properties 
are given as £=6.825x10° and v=0.3. The shell is subjected to concentrated 
forces as shown in Fig. 7. Only 1/4 of the shell is modeled and symmetry 
bouncary coaditions are imposed on two sides. In addition the * node is 
constrained in the Z direction. Convergence 1s studied by varying the mesh 
sive, Input and outut are given for an 8x@ mesh. A theoretical lower bourid 
for the radiat deflection at the center in the case where the hole at the 
center is not present rs .0924. A value of .0940 was used for normalization 


of tne results. 
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INPUT TS PROVIDED 


PROBLEM #8 
REGULAR SYMMETRIC CROSS-PLY LAMINAE 


FIGURE & 


A composite square plate of jeagth 5 modeled as a symmetr’: cross-ply laminate 
[0, 90, 0) with four elements ts shown in Fig, &. The thickness of each ply 
is ,02222, The miteria] properties are given as E,=2.0xlu’, E25 0x10, 

Vv, 520.25 and 6,,*2.5xiG" and are the same for eech ply. Three“corners of the 
pt Ste are pinneg and a unit force in the -7 direction ts applied at the free 
corner to simulate © pure twise loadiag. lhe theoretical selution for the Z 
deflection at the free corner is -3.750x10°~, Theoretical resuits are also 
given for T for tie element cantrining the free corner in Table 3, 


1 1.50.0 
? 0.0 
3 50.0 | 
TABLE 2 


INPUT YS PROVIDED 


PROBLEM #¥ 
REGULAR SYMMETRIC CROSS-PLY LAMINATE 
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FIGURE 9 


A quarter model of a composite square plate modeled as a symmetric cross-ply 
laminate is shown in Fig. 9. The lTength of each side 15 1.0, and the thick-. 
ness of egch piv is .000666. The material properties are given as E,=2.0x10', 
yeh 0x107, “70.25 and G,.=2.5x10" and are the same for each ply. “The full 
ofate is simp! § supperted dha is subjected to a uniform pressure Joad of 
-1.0x10°°. The finite elemcnt model contains 25 nodes and 16 elements, The 
theoretical solution for the 2 deflection at Grid 25 (center of the plate) is 
-1.836x19 ~, Theoretical resuits are also given for the stresses in Jayers 1 


and 3 in element 17 in Table 4. 
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TABLE 4 


INPUT AND OUTPUT ARE PROVIDED 


PROBLEM #10 


COMPOSITE OPEN TUBE 


A finite element model of a composite open tube modeled with the symmetric 
lay-up £45, -45, 0, 90, 90, 0, -45, 45] is Shown in Fig. 10. The model 
contains {44 nodes and 128 elements. The length of the tube is 80, the 

radius 15 50, and the thickness of egch layer is .24,. The material properties 
are civen as £,273.8x10°, £,=3.75x10°, V,,=0.4, and Gy=1.74x10"- The tube 

is Subjected t6 a uniform pfessure Toad be 10.5. The’fheoretical solution 

for the hoop loading, FY, for element 8 is 525. Results from MSC/NASTRAN 

are given in Table 5 for the layer stresses in element 8. 


i a 
1 2.519E2 | 1.74061 2.27361 
2 2.49462 | 1.74861 | -2.273E1 
3 -?.25762 | 3.2886] | } .L08E-2 
4 7 PbOE2 b 2cGbIEG: 4 2. eT EA? 
6 196062 | 2.65760 F 5.551E~2 
6 -2 25502 | 3.2316) | -2.a8i£-2 
] 953562 | 1.74161 | ~2.274E1 
i | @ 47862 | 1.7621 | 2 37AE] 
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INPUT IS PROVIDED 


PROBLEM #11 
STRAIGHT BEAM TEST 


FIGURE 11 


A finite element model of a cantilevered beam modeled as a laminate configura- 
tion [0, 0, 0, OJ in a simulation of equivalent isotropic properties is shown 
in Fig. 11. The length of the beam is 6, and the thickness of each layer is 
-25. The material properties are given as £=10.0x10° and v=.3. The beam is 
subjected tu two loading conditions as follows: 


CASE 1: A unit force is applied at the tip in the X direction. (Extension) 
CASE @: A unit force is applied at the tip in the Z direction. 

For Case 1 the theoretical solution for the X deflection at Grids 13 and 14 is 
3.0x10° >. For Case 2 the theoretical solution for the Z deflection at Grids 


13 and 14 is .4320. Theoretical results are given for the bending moment 
distribution from the free end to the cantilevered end in Table 6, 


TABLE 6 


For Case 2 results from MSC/NASTRAN are given in Table 7 far the direct layer 
bending stress in element 6. 


1.23854 

4.12653 
-4.126E3 
-1.23854 


TABLE 7 


INPUT AND OUTPUT ARE PROVIDED 


PROBLEM #12 


COMPOSITE SHELL ROOF 


FIGURE 12 


A finite etement model ot a composite shell rocf modeled with the symmetric 
jayeup [45, -45, 15, -15, -15, 15, -45, 45] is shown in Fig. 12. The length 
acd radius of the shell are 25, and tne thickness of each jayer is .03125. 
The material properties are given as E,=2.0x10°, £,70.5x10°, V,,=0.25, 
G)p=2.5x10°, and G,o=G,,=2.5x10". Thelshell is subjected to a Gniform 
ph&ssure of 90.0. Reshfts from MSC/NASTRAN are given in Table 8 for the 
radial deflection at selected nodes. 


of 
34.) ~1 0662 


35 ~1.3441 
| 36 -1.6074 
43 ~1.3267 
44 ~1.6739 
45 ~Z 6079 


Aa ehh ee 


TABLE 6 
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PROBLEM #13 
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Oo, 0, 90, 
elements. 
er each layer 
E,=1.6x108, 
plate are conside 
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The length of the plate is 20, 
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COMPOSITE RECTANGULAR PLATE 
132 
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FIGURE 13 


A composite rectangular plate mcdeled with the symmetric lay-up [45, 
contains 
the width is 
proverties are given as 


] is shown in Fig. 13 The model 


The material 
.65x10° and G,,=G,,=4.0x16°. 


.0236. 
Gy 


red. 
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dui the -X direction at the %-15.0 ede 
Cy 2 The plate as modeled with bending pro 
unifGrm pressure icad of 2 in the -Zd 

Gi 


Hi 


-45, 90, 
36 nodes and 25 
and the thickness 
Melis. Selo, 

of the 


beer 


a 
i 


loa 


ammonia EN AR a Mier th a MSTA RC NS RCA NY NN in, am i eae il fai Ce UO SS OR Re ecenes” SERS ei aCe Rao OT es “IE RR oie ee 


CASE 3: The plate is modeled in membrane and bending and ts subject to the 
loads defined in Case 1 and Case ?. 


CASES 1 through 3 are defined for a static analysis, RF #1. 


CASE 4: Normal modes analysis using RF #3. The first ten natural frequencies 
are found, 


PROBLEM #14 


AY COMPOSITc SANDWICH PLATE 


FIGURE 14 


A sanawich plate with an isotropic core and c.nposite skins modeled with 
the symmetric lay-up [45, -45, 90, 0, 0, 90, -45, 45] is shown in Fig 14. 
The model contains 35 nodes and 25 elements. The length of the plate is 20, 
the width is 15, the thickness of each layer is .0235, and the thickness of 
the core is 1. The material properties of the composite skin are u en as 
E,718.5x10", tp21.Gx10°, v0.25 and 6)o°.65x10°, 6, 7=G.7=4.0x10°. The 
material properties of the core are given as G=3.0x10'. The plate 15 


subjected f2 a uniform pressure load of 2. 
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AFPEPDIX C: CASE CONTROL AND BULK DATA CARDS 
PERTAINING TO THE QUAD4 ELEMEWi 


Case Control Data Card ELSTRESS 


Description: 


ML NTEN OA a ERE SEY SMEVIAG) INA a pe * A ee 
oO 8S a em = — 


« Limwent Stress Output Request 


Requests farm and type of element stress output. 


Format and Exampie(s): 


ELSTRESS 


ELSTRESS » 
ELSTRESS «# 


ELSTRESS(SPRT2, PRINT, PUNCH, PHASE) = 25 
Option Meaning 


SORT 


SORT2 


PRINT 
UNCH 


EXTREME or 
LAYER 


REAL or 
IMAG 
PHASE 
ALL 

n 


NPNE 


Remarks: 


EXTREME 


SBRTI 
UAVER™ =, TAG 


PRINT 
« BUNCH , 


ALL 
od h 
ASPRINT NBNE 


ALL 


Output will be presented as a tabular listing of elements for each Toad, 


frequency, eigenvalue, or time, depending on the rigid format. S@RT1 is not 


available in Transient problems (where the default 4% S@RT2). 


Output will be presented as a tabular listing of load, frequency, or time for 
each element type. SPRT2 is available only in Static Analysis, Transient end 
Frequency Response problems. 


The printer will be the output media. 
The card punch will be the outsut media. 


Requests stresses to be calculated at the extreme (top and bottom) fibers of 
8 plate element or, for composites, the stresses for each layer. (See Remark 8) 


Requests real and imaginary output on Complex Figenvatue or Frequency Response 
problems. 


Requests magnitude and phase (0.0° « phase < 360.0°) on Complex Eigenvalue or 
Frequency Response problems. 


Stresses for all elements wil? be output. 


Set identification of a previously appearing SET card (Integer > OG). Only 
Stresses for elements whose identification numbers appear on this SET card will | 


be vutput. 


Stresses for no elements will be output. 


Both PRINT and PUNCH may be requested. 


ALL in Trans:ent and Frequency response problems geners! ty 


Ax output request fr- 
produces Yarge amounts of printout. An alternative to this would be to define a SET 
| 
\ 


of interest. 


In Stetic Analysis or Frequency Response problems, any request for SORT2 cutout a 
causes All output to be SPRT2. 
{ 

I 


oo, 


CSTRESS is an alternate farm and is entirely equivalent to STRESS. 


r 


ELSTRESS © NBRE ailows overriding an overall request. 


1 
(Continued) 


2.3-18 (8; 1Oo/8? } 


i1& 
S 


ELSTRESS (Continued) 


‘ sd 6. if element stresses in material coordinate system are desired (only for TRIAI, 
TRIA2, QUAD] and QUAD2 elements and only in Rigid Format 1), the parameter STRESS 
(see the description of the PARAM bulk data card fn Section 2.4.2) should be set to 
be e positive integer. If, in addition to element stresses {mn material coordinate 
system, stresses at the connected gric points are also desired, the parameter STRESS 


should b2 set to 6. 


zs @ 7. When LAYER fs selected, inaividual layer stirvsses and/or failure indices will be 
output. 
8. The option £XTREME and LAYER is gniy applicable for the QUADS wlesment. 


£.3-1Ba (BS ECrEY,; 


ens hte nee oe UnenenaSas th -tedlmi ue. \ kteleanae “ SAREE ALPS Net CU IONS NRK 
AA, SOA! Sa Ue ek Rael nelN fe SiemvedmaRK sane i ‘ARERR CA orci, salen: 
a x é 


bescription: 


STRESS 


Option 


EXTREME or 
LAYER 


REAL or 
IMAG 
PHASE 
ALL 

n 


NONE 


Remarks: 


STRESS = ALL 
SYRESS(S@RTL, PRINT, PURCK, PHASE) = 15 


tab 


Cese Contro’ Bata Card STRESS - Element Stress Output Request 
Requests form and type of element stress output. 


Format and Exemple(s)}: 
sort: PRINT —-ERTREME = REAL ALL 


SORTC =, BUNCH , LAYER” , VRAG Jie 4 on 
NOPRINT PHASE NONE 


Peaning 


Output will be presented as a tabular Tisting of elements for each toed, 
frequency, eigenvalue, or time, depending on the rigid format. SBRT1 is not 
available im Transtent problems (where the default is S@RT2), 


UYutput will be presented a3 a tabular isting of load, freouency, or time for 
each element type, SORT? ts available only in Static Analysis, Transient and 
Frequancy Response problems. 


The printer will be the output mecie. 
The card munch wilt be the output medta. 


Requests stresses to be calculated at the extreme (too and bottom! fibers of 
a plate element or, for composites, the stresses for each layer. (See Remark 8) 


Requests real and imaginary output cn Complex Eigenvaiue ov Proaequancy Response 
problems. 


Requests magnitude and phase (0.0° s phase « 360.0°) on Comples tfgenvatue or 
Frequancy Response problems. 


Stresses for all elements will be output. 


Set identification of a previously appearing SET card (integer >» 0). Only 
stresses for eiements whose identification aumbers appear an this SFT card will 
be outpet 


Siresses for no elements wil! be output. 


Both PRINT ang PUNCH may be requested, 


An output request for ALL on Tranafent and Frequency response problems generally 
produces Varge amounts of priuntaut. 
of interest, 


In Static Analysis or Frequency Response problems, any request flow MERTY cetout 


causes Stl outout to be SBRT2. 
STRESS is an alternate form end fs entirely squdwalent to EGSTRESS 
STHESS « NENE allows overricisg an overai) request. 


{Corntiaued) 
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SURESS (Continued) 


6. If element stresses in materval coordinate system are destred {only for TRIAI, 
TRIAZ, QUAD] and QUADS elements and oniy fm Rigid Format 1), the parameter STRESS 
{see the description of the PARAM bulk data card in Section 27.4.2} should be set te 
be a positive integer. If, 4 edition to elament stresses {fn materte] coordinate 
system, stresses at the connected grid points are also dattred, the garameter STRESS 
should be set to 0. 


7, When LAYER is selected, findividunl fayer stresses and/or Failure dadieas will be 
Output. 


B. The option EXTREME and LAYER is only applicable for the QUAD4 element. 


Input Data Card CQUAD4 Quadrilateral Element Conrection 


Sescription: Defines a quadrilaterial plate element (QUAD?) of the structural model. This 15 an 
fsonarametric ee at element, with variable element thickness, Jayered 


composite material and thermal analysis capabilities. 


Forwut end Exumole: 


Field Contents 
£10 Element identification number (Integer > 0) 
PID Icentificetion number of a PSHELL entry (Defaule is EIR; (Integer » 0) 
For composites, see Remark 5. 
G3 Srid point identification numbers of connection points (lateger > 0) 
20 Offse: uf the element reference plane from tie plane of gria potnts (Real or 


blank, see Remark 3 for default) 


Th Materts) property orfentation specification (Real er blank; or O s Integer 
« 1,000,000). [f Real or blank, specifies the material property orientation 
angle in degrees. If Integer, the orientation of the material x-axts is 
along the projection onte the chane of the alement of the x-axfs of the 
coordinate system specified by the integer valve. 


Mombrane thickrass of element at arid potets Gi (Peal cr blank, see Remark 


4 vor default). 


Remarks: 1. Tne QUADS geometry, coordinate systems and numbering are shown in the figure below: 


ay 


Ny ; oss mene 


2. Element identificet.on wumbers must be unfque with respect to all other element 


identiffearton numbers . 
{Continued} 
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CQUAD4 (Continued) 


3. The mterial coordinate system (TM) and “he offset (270) may also be provided on the 
PSHELL entry. The PSHELL data will be used if the corresponding field on the CQUAD4 
entry 1s blank. 


4. The T4 are optional, if not supplied they will be set te the value of T specified on 
the PSHELL entry. In such cases, the continuation entry ts not required. 


5. For composites, a PCOMP, PCOMP], PCOMP2 card can be used instead of a PSHELL card. 


2.4-87b (&/10/67) 


BULK DATA DECK 


Input Data Card MAT2 Material Property Definition 


Description: Defines the material properties for linear, temperature-independent, anisotropic 


matertals. 


Format and Example: 


7 2 3 4 g 6 7 8 ] 10 


;0.956 | Asc 


+def_ MCSID | : | 
+EF 1008 me ose = -_ 


Field Contents 

MID Material identification number (Integer > 0) 

Gij The material property matrix (Real) 

RH Mass density (Real) 

Ai Thermal expansion coefficient vector (Real) 

TO Thermal expansion reference temperature (Real) 

GE Structural element damping coefficient (Real) 

ST, 5C€, SS Stress limits for tension, compression and shear (Real) (Used only to compute 


margins of safety {n certain elements; they have no effect on the computational 


procedures ) 
MCSID Material coordinate system identification number (integer > 0 or blank) 


Remarks: 1. The material fdentification numbers must be unique for all MATT, MAT2 and MAT3 cards. 


2. MAT2 materials may be made temperature dependent by use of the MATT2 card. 


X. The mass density, RH®, will be used to automatically compute mass for al} structural 
elements except the two-dimensional bending only elements TRBSC, TRPLT and QOFLT. 


4. The convention for the Gy in fields 3 through B is represented by the matrix 
relationship. 


42) ae e 
&, 


5S. MCSIN (> 0) is required if stresses or strains/curvatures are to be computed in 


material coordinate system. This is applicable only for TRIAT, TRIA?2, QUAD), and 


QUADE elements. 


2.4- 169 (12/23/78) 


Input Data Card MATS Orthotrepic Plate Material Property Definition 


Description: 


Format and Example: 


Defines the material property for an orthotropic material for plate elements. 


Fle 


Contents 


Material identification number (Integer » 0) 

Modulus of elasticity in the material x and y directions (Real #¢ 0.0) 
Poisson's Ratio (Real) (See Remark 5) 

Linear In-nlane shear modulus (Real > 0.0) 

Transverse shear modulus for shear in X-Z plane (Real) 

Transverse shear modulus for shear in Y-Z plane (Reai) 

Mass density (Real) 


Thermal expansfon coefficfents in the material x and y directions 
(T, Real > 0,0) 


Thermal expensson reference temperature (AC, Reu)) 

Allowable stresses/strains in tension and compression, respectively, 4m tne 
material x direction. Required if faflure index calculation is desired. 
{XT, Real > 0.0) (XC, Real) (Default value for XC is XT) (See Remark 3) 
Allowable stresses/strains in tension and compression, respectively, fn the 
material y direction. Required if failure index calculation {s desired. 
(YT, Real > 0.0) (YC, Real) (Default value for YC 1s YT) (See Remark 3} 
Allowable stress/strain for in-plane shear (Real > 0.0) (See demark 3) 
Structural damping coefficient (Real) 


Tsai-Wu interaction term (Real) (See Remark 4) 


(Continued) 
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Remarks: 1. 


ba 


uw 


MATE (Continued) 
Material coordinate systems are defined by the plate element connection entries on 
the CQUADS card, 
The stress-strain relationship defined by this data is: 
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Fields XT, AC, YT. YC and S ere used only for composite materials when felture 
calculations are requested with PCOMP, PCOMP) or PCOMP2 Bulk Date entries. 
Ailowables represent stresses except when the maximum stvaia failure theory {s used, 


The Fl2 fieid is used only for composite materials when the Tssi-Wu fadlure theory 
{s used and failure calculations are requested. 


NU12 fs Pottson’s Ratio (e,/ey for unfaxtal Yoeding in l-directiorn). fote that 
NUZI tEo/ fy unjaxial Toading in 2-direction, {s related to NUI2, El and E2 ty 
the relationship, (KUI2) (£2) =» (NU2T) (ET). 
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Input Data Card PCOMP Layered Composite Element Property 
Description: Oetines the properties of a n-ply laminated composite material. 


Format and Example: 


COMP PID 


eee err aS wenn ENRON 1 


COMP 106 


anak Nee 


THI T sours ETC, er 
45.0 w.oR 
ae = Seeee ee See Rae 


1a Contents 

PID Property identification number (1,070,900 » Integer > 0) 

AG Offset of the element reference planw from the plane of arid points (Real 
or blank, see Remork 2) 

NSM Non-structural mass per unit area, (Real) 

SBOND Allowable shear stress of the Lending material. (Real » 0.0 or blank) 
Required if fatlure theory is used. 

rT Failure theory, one of the strfngs “NILL", “HOFE", “TSAIS, "STRESS", or 
"STRAIN". See Remark 4. (BCD or blank} 

LOFT Lamination generation option, cne of the strings, "ALL™, "SYM", SME", or 
"SYMMEM' Sea Remark 3. (BCD or blank) Defeult ts all. 

MIDS Material identification number of the en Yayer. (Integer > 0 or blank) 

TH Thickness of the ght Teyer (Real, > 0.0 or Blank} 

Thi Angle between the longitedinal direction of the fibers of the on leyer 


and the material X-axis. (Real or b)ank) 

soumd Stres* output request for ith later, one of tre strings "YES" or "NO". 
(Pefault is “NO") 

Renarks: 1}. Yhe plies are numbered from to tuo rp beginning with the bottom layer. 


2. The offset (20) ts used only when the corresponding field on the CCRAD4 Bulk Dets 
entry referencing this property are blank. 


{Continued ) 
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PCOMP (Continued) 


SBOND 4s required {f bonding material failure index calculations are desired. 


&. The failure thecry if used to determine the element failure on & ply-by-ply 
besis. The available theories wre: 


HILL - Hill Theory 

HOrF « Hoffman Theory 

TSAI - Tsai-dy Theory { 
STRESS ~ For Maximum Stress Theory 


STRAIN ~ For Maximum Strain Theory 


5. To minimize input requirements several lamination options (LIFT) are avetiable. Al 
4ndicates that every ply fs specified, SYM indicates that ply leyup ¢5 symratric 
about the center ply and that the plies on one side of the center Hine ere 
specified  SYMMEM indicates @ symmetric layup cf membrane only plies. 

o. The matertal properties, MID{, may reference only MAT], MAT2 and HATS Bulk Date 
entries. 

7 Yt any of the MIDI, Ti or THI are blank, then the Tasy roneblank values specified 
for #ach will be used to define the values for the ply 
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Input Data Card PCOMP] Layered Composite Element Property 


Description: Defines the properties of a n-ply laminated composite material where all plies 
are composed of the same material and are of equal thickness. 


Format_and Example: 


pT eT Toe eo Der 
poor | too j 0.5] 17 {500 f swam] 200 | 026 | sm | ac 


Field Contents 

PID Property identification number (1,000,000 > Integer > 0) 

20 Offset of the element reference plane from the plane of grid points (Real 
or blank, see Remark 2) 

NSM Non-structural mass per unit area. (Real) 

SBOND Allowable shear stress of the bonding material. (Real > 0.0) 

FT Failure theory, one of the strings "HILL", “HOFF”, "TSAI", “STRESS”, or 
"STRAIN". See Remark 4. 

MID Material identification number for al] layers. (Integer > 0) 

LOPT Lamination generation option, one of the strings, "ALL", "SYM", "MEM", or 
"SYMMEM". See Remark 5. 

TPLY Thickness of all layers (Real, > 0.0 or blank) 

THI Angle between the longitudinal direction of the fibers of the jth layer 


and the material X-axis. (Real or blank) 


Remarks: 1. The plies are numbered from 1 to n beginning with the bottom layer. 


2. The offset (20) is used only when the corresponding field on the CQUAD4 Bulk Data 
entry referencing this property is blank. 


3. SBOND fs required if bonding material failure index calculations are desired. 


4. The failure theory fs used to determine the element faflure on a ply-by-ply 
basis. The available theories are: 


HILL - Hill Theory 

HOFF - Koffagn Theory 

TSAI - Tsai-Wu Theory 

STRESS - For Maximum Stress Theory 
STRAIN - For Maximum Strain Theory 


(Continued) 
2.4-223c (8/10/87) 
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PCOMP) (Continued) 


To miniinize fnput requiremerts several Taminetion options (LOPT) ave available. Atl 
{Indicates that every ply ts specified, SYM indicates that ply layup {s syametric 
about the center ply and that the plies on one side of the center line are 
Specified. SYMMEM indicates a symmetric layup of membrane only plies. 


The material property, MID, may reference only HAT], MATZ and MAT& Bulk Data 
entries. 


Input Data Card PCOMP2 Layered Composite Element Property 


Description: Defines the properties of a n-ply laminated composite material where all plies 
are composed of the same material. 


#ormat_and Example: 
1 2 3 4 § 


6 7 
SBOND FT =: «IO 
6.43 TSAI 


Field Contents 

PID Property identification number (1,000,000 > Integer > 0) 

£0 Offset of the element reference plane from the plane of grid points (Keal 
or blank, see Remark 2) 

NSM Non-structural mass per unit area. (Real) 

SB0ND Allowable shear stress of the bonding material. (Real > 0.0) 

4T Failure theory, one of the strings "HILL", "HOFF", "TSAI“, "STRESS", or 
"STRAIN". See Remark 4. 

MID “Material identification number for all layers. (Integer > 0 or blank) 

LOPT Lamination generation option, one of the strings, "ALL", "SYM", "MEM", or 
"SYMMEM". See Remark 5. 

Ti Thickness of the qth Tayer (Real, > U.0 or blank) 

THi Angle between the longitudinal] direction of the fibers of the jth layer 


and the material X-axis. (Real or blank) 


Remarks: 1. The plies are numbered from 1 to n beginning with the bottom layer. 


2. The offset (Z0) is used only when the corresponding field on the CQUAD4 Bulk Data 
entry referencing this property is blank. 


3. SBOND is required if bonding material failure index calculations are desired. 


4, The failure theory is used to determine the element failure on a ply-by-ply 
basis. The available theories are: 


HILL - Hill] Theory 

HOFF - Hoffman Theory 

TSAI - Tsai-Wu Theory 

STRESS - For Maximum Stress Theory 
STRAIN = For Maximum Strain Theory 


(Continued) 
2.4-223e (8/10/87) 


129 


@ 
PCOMP2 (Continued) 
§. To minimize input requirements several lamination options {LOPT) are available. ATT € 
indicates that every ply is specified, SYK indicates that ply layup 13 symmetric 
about the cenrer ply and that the plies on one side of the center line are 
epecified. SYMMEM indicates o symmetric layup of membrane only plies. 
6. The material property, &I9, may reference only MATI, MAT2 and MATS Bulk Data 
entrigs. 
7. If any of the Ti or TH1 are blank, then the last non-blank values specified for each a 
wilt be usec to define the values for the ply. 
@ 
4@ 
é 
@ 
| 
@ 
@ 
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aInput Data Card PLOADS Fressure Loeds on Face of Structural Elements 


fescription: Defines @ load on a face cf & QUADS elewent. 


” 
e 


Format and Example: 


awa NNN 


+h4 Cio 

“i oe fe aie 2375 

Fleld Contents 

Sec Load set {dentification number (Integer > 0) 

EID, EL,E2 Etement idencification number (Integer > 0, El « tj 

2) Pressure at the grid point defining the element face (Real or blank) 
elo Coordinate system identification number (Integer » 6) 

Ai Components of a vector in system CiD that defines the direction of the 


prassure (Real; 


Remarks: 1. For the plate element QUADS, 1¥ the continuation entry fs not given, the direction 
of the pressure is norms! to the element in the element Z direction, If only Pi is 
given, the pressure {5 assumed to be uniform over the element surface. 


?. If the loaded surface of an element is curved. and a direction vector 1s not 
specified, the direction of the pressure mey vary ever the surface. The pressure 
intensity 18 the load per unit surface area. 


3. Equivalent grid point Tuads ere computed. A unfform pressure teed nor reselt in 
@qua!l grid point loads. 


78-2395 (8/10/87) 


at 


Hulk Data Entry PSHELL Shell Element Property 


cr’ Defines the membrane, bending, transverse shear, and coupling propert-.es of 
the “YAD4 digell element. 


Formt and Example: 


i 2 3 
peers pemegreece 
PSHELL & PID M101 


PSHELL 203 204 


ae Saree 


Fieid Contents 

PID Property ¥dentiticatign number (Integer > 0) 

Mig) Material identification number for membrane (Integer > 8 of dlant) 

a Default value for membrane thickness (Real » 0.9) 

M102 Material identificetion number for bending (Intever » OG or blank) 

\z/T Bendlng stiffness parameter (Real or blank, default = 1.2) 

HID3 Material identification number for transverse shear (Integer » 0 or biank, 
must be blank unless MID2 > 0) 

TS/T Transverse shear tiickness divided by membrane thickness (Real or Biank, 
default = .833333) 

NSM Nonstructural mass per unit area fkeal} 

71,22 Fiber distances for stress computation. The arsitive direction: fs determined 


by the righthand rule and the ovder fn which the grid potiuts ere listed on the 
(Real or tlank, defaults ave -T/? for Z] end 43/2 fer Ze). 


connection entry. 


MID4 Material identification number for membrane-bending coupling (Integer >» 0 or 
otank, must be blark unless MiDL > u and MID2 » 3, may not squall KIDI or M102) 


MCSTO Identificatfon number of matertel coordinate system (Real or lank, or 
(integer & 0) (Sec Remark 32) 


SCSI0 Tientification wumber of stress coordinate system (Real or blank, or 
{Integer z G) (see Nemes 11) 


rad) Offset of the element reference plane “rom the prane of grid points. 


aor Stank, defzytt = 0.0) (See Remark a2) 


(Continued) 
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PSHELL (Continued) 


Memarks: 12, All PSHELL property entries gust have unique identification numbers. 


2. The structural mass fs computed from the density using the merbrane thickness and 
wenbrane material properties. 


3. The results of leaving any HID field blank are: 


MID) No membrane or coupling stifiness 

MIDZ Ho bending, coupling, or transverse shear stiffresy 
HIDS No transverse shear flexibility 

MIDée Wo membrane-bending coup! tag 


4. The continuation entry is hot. required. 
5. Strvctuval damping. when needed, is obtained from the MIDI seaterca. 


6. The WI04 field should be laft blank if the materdal properties are symmetric with 
the middle surface of the shelT. 


7. For struatural problems, PSHELL entries may reference MATL, MAT2 or WATB material 
prone date, 


8B. I% the transverse shear material, MID3, references MAT? date, then 633 gust be zero. 
If AIO2 references MATS data, then G1,2 and G2,Z must not be zero. 


9. For Seat transfer problems PSHELL entries may reference 4AT4 or MATS material 
property data. 


1U. Jf MCSID/SCSID ts left blank (0.0) or ts real, it is considered to be the angle of 
rotation of the X axis of the materfal/stress coordinate system with respect to the 
Rh axis of the element coordinate system in the XY plane of the latter. If Integer, 
the orientation of che matertal/strecs x-axis 1s along the projection of the x-axis 
of the specified coordinate system onts the x-y plane of the element system. The 
value of MCSID is the default value for the 1% field on the CQUAD4 Bulk Date 
entries. 


Ll. The value of 20 is the defsult value for the corresponding field on the CQUAD4 Bulk 
Nata entries, 
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APPENDIX D: WASTRAN INPUT DATA FOR SAMPLE 
PROBLEMS 1-6, 8, 9, 11, 13c, 14 
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PROBLEM #9 
REGULAR SYMMETRIC CROSS-PLY LAMINATE 
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FIGURE 9 


A guarter model of a composite square plate modeled as symmetric cross-ply 
Tem siete ts snown in Fig. 9. The length of each side is 1.0, and the thick, 
ness of each ply is .C90666. The material properties are given as E,=2.0x10 ; 
Eoe5 0x10", % 80.25 and G,462.5x10° and are the same for each ply. “The full 
plate is,simp!§ supported Aa 1s Subiected to a uniform pressure load of 
-1.0x10 ". The finite element model contains 25 nodes and 16 elements. The 
treoret{ical soiution for the Z deflection at Grid 25 (center of the plate) is 
~1.836x10°°, Theoretical results are also given for the stresses in layers j 
and 3 in element. 17 in Table 4, 
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PROBLEM #12 


COMPOSITE SHELL ROOF 


FIGURE 12 


A finite element model of a composite shel} 
lay-up (45, -45, 15, -15, -15, 15, -45, 45] 


roof modeled with the symmetric 
is shown in Fig. 12. The length 


and radius of the shell are 25, and the thickness of each jayer is .03125. 


The material properties are given as E,=2.0x10°, & 


G,.=2.5xl0 , and G, = =2.5x10°. Theshell 
ptassure of 90.0. !Resatts from MSC/NASTRAN 


radial deflection at selected nodes. 


=0.5x10 7) V1 250.25. 
is subjected to a uniform 
are given in Table 8 far the 
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PRORLSM #15 


COMPOSITE RECTANGULAR PLATE - UNSYMM:TRIC CROSE-PLY LAMINATE 


UNSYMMETRIC FIGURE Al 
CROSS-PLY LAMINATE 


One half cf a composite rectangular plate modeled with the unsyametric crose- 
ply layup (0°, 90°) is shown in Fig. Al, The plate is mcucled with 15 nodes 
aad 8 elenente. The length of the plate is 5", the width is 1", and che 
thickness of each layer is .005". The material properties matrix on the 

MAT2 card te given as 6), = 1.94x10’, Gy, = 4.66x10", Gy, = 1.33x10° end 
G33 * 8.319. The thermal expansion coefficient vector is defined by 

Ay = ~Exto7® Ay = t8x107° and A,»> = 0. A temperature field was detinec 
for the plate by an average temperature of 230° over the cross-section of 
each element on the TEMPP] card. The input deta, the calculated thermal 


luedirg, snd the dispiacements due to the viermal loadi 2 are also provided. 
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